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Princeton  VLSI  Project:  Semi-Annual  Report 

R.  J.  Lip  ton 


1.  Introduction 

We  have  been  officially  underway  at  Princeton  tonce  August;  hence,  this 
report  actually  covers  about  three  months  of  activity Si There  are  three  major 
aspects  to  our  project:  ALI,  Census,  and  Test  ing 


; 


2.  AL1 


<~>  ALI.  our  procedural  language  for  VLSI  design  and  layout,  is  now  up  and  run- 
ning,{4^>  It  has  already  been  used  by  Dobkin  and  Drysdale  to  redesign  a  divider 
they  previously  designed  on  a  graphics  system  at  Xerox.  In  addition.  LaPaugh's 
VLSI  class  has  already  used  ALI  for  their  VLSI  projects.  Dobkin  is  currently 
beginning  to  explore  ways  to  create  graphics  interfaces  to  the  ALJ  system. 

ALI2  the  second  version  of  our  layout  system  is  coming  along  well.  Valdes 
has  already  fully  defined  the  new  language  and  implementation  is  now  unker  way 
[7]  AU2  differs  from  ALJ  in  two  essential  ways.  First,  it  is  based  on,  iwhat  we 
believe  to  be  a  much  cleaner  set  of  primitive  and  constructs.  For  instance,  it 
does  not  have  the  "shuffle  property"  so  many  design  languages  do-  by  this  we 
mean  that  the  order  of  placement  commands  does  matter.  Secoryl  AL12  gen¬ 
erates  far  fewer  constraints  than  ALJ  does.  This  is  of  course  critical  if  such  con¬ 
straint  based  systems  are  to  be  able  to  handle  large  complex  layotits 

Plans  are  already  underway  on  how  to  best  exploit  the  unique  features  of 
AL12.  Ramachandran  has  just  finished  a  study  of  the  cost  of  increasing  the  size 
of  drivers  in  order  to  speed  up  circuits  [5].  She  assumes  that  sizes  of  transis¬ 
tors  can  be  changed  but  that  the  layout  cannot  be  restructured.  Under  these 
assumptions  she  can  tightly  bound  the  worst  case  cost  in  terms  of  area  of  mak¬ 
ing  the  delays  on  all  wires  a  "constant".  Since  A1J2  allows  a  designer  to  easily 
change  the  sizes  of  wires  and  transistors  we  feel  that  such  results  are  important 
Already  it  is  common  place  for  designers  to  size  transistors  for  speedup  in  an  ad 
hoc  way  [6] 

Also  Vijayan  and  Wigderson  have  isolated  a  number  of  new  layout  problems 
that  arise  naturally  when  one  considers  the  implementation  of  AL12  [8.9].  Ail 
these  problems  concern  the  embedding  of  "rectilinear"  graphs,  these  are 
graphs  where  each  edge  is  connected  to  the  "left”,  "top",  "right", or  "bottom"  of 
each  vertex  They  are  currently  exploring  the  computational  complexity  of  a 
variety  of  layout  problems  here.  For  example,  they  can  quickly  recognize  those 
graphs  that  have  planar  embeddings;  moreover,  they  can  also  quickly  find  such 
embeddings  when  they  exist. 
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Census 

'A  The  census  language  is  a  new  way  to  express  parallel  algorithms  that  use  a 
fairly  loosely  coupled  method  of  control,  Work  is  under  way  to  understand 
the  limits  and  powers  of  such  languages.  North  is  beginning  to  identify  those 
problem  areas  that  can  be  successfully  mapped  onto  the  census  language.  We 
are  also  thinking  about  implementations  of  census,  but  no  implementation  is  yet 
started.^ 

Census  has  been  looked  at  from  a  theory  point  of  view  by  Chandra,  Fortune, 
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and  Lipton  [2].  They  have  been  able  to  get  tight  upper  and  lower  bounds  on  the 
size  of  boolean  circuits  with  unbounded  fan-in.  Unbounded  fan-in  circuits  not 
only  model  an  important  class  of  census  computations  but  they  also  model  cir¬ 
cuits  such  as  PLA's.  For  example,  it  is  possible  to  construct  a  circuit  that  adds 
two  n-bit  numbers  in  constant  time  and  whose  size  is  approximately  linear  in  n. 
We  are  beginning  to  examine  the  feasibility  of  using  these  ideas  in  real  VLSI 
designs. 

Testing 

Our  work  on  testing  divides  into  two  areas.  Arden  is  beginning  to  work  on 
first  silicon  testing  especially  with  respect  to  scanning  electron  microscopes 
(SEM).  He  plans  to  be  on  leave  this  spring  and  work  with  the  SEM  group  in  Mun¬ 
ich.  In  the  future  we  expect  that  we  will  be  able  to  use  a  SEM  that  the  Princeton 
Siemens  group  is  about  to  get. 

LaPaugh  and  lipton  have  begun  to  work  in  the  area  of  production  testing 
They  have  already  successfully  been  able  to  completely  characterize  the  testa¬ 
bility  of  "prefix  computations".  Prefix  computations  arise  naturally  in  a  number 
of  places,  for  instance,  in  carry-look-ahead  adders,{4^>  Such  characterizations 
link  the  self-test  of  these  computations  with  classic  semi-group  theory. 

Also  work  is  under  way  on  a  new  self-test  strategy  which  we  call  "toggle 
search".  This  method  first  generates  the  vector  of  all  0‘s;  next  it  randomly 
changes  one  of  the  0’s  to  a  1;  it  repeats  this  until  all  0's  have  been  toggled  to  l’s, 
then  the  entire  process  begins  again.  Toggle  search  and  its  generalizations  are 
well  suited  to  many  test  environments  where  one  bit  of  an  input  can  be  changed 
more  quickly  than  a  whole  input  vector.  Already  we  have  empirical  evidence  of 
the  superiority  of  toggle  search  over  other  methods.  Colleagues  at  IBM  Watson 
Research  have  tested  toggle  on  examples  of  about  one  thousand  gates  and  found 
that  it  is  about  three  times  faster  than  standard  methods.  For  example,  it  took 
about  600  thousand  vectors  versus  1.8  million  to  achieve  100%  fault  coverage  on 
one  piece  of  random  control  logic.  We  plan  further  experiments  to  further  vali¬ 
date  these  results. 

Finally,  we  have  also  found  a  way  to  transform  any  combinational  logic  cir¬ 
cuit  into  one  that  is  easy  to  test.  Here  by  easy  to  test  we  mean  that  we  can 
detect  a  very  large  class  of  physical  faults.  The  penalty  for  this  transformation 
is  that  the  number  of  gates  can  increase  by  as  much  as  100%.  We  plan  this  com¬ 
ing  year  to  carefully  explore  this  new  method.  In  particular,  we  wish  to  both 
understand  the  cost  of  the  method  and  the  class  of  faults  it  can  and  cannot 
detect. 
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Abstract:  The  first  component  of  a  VLSI  design  environment  being  built  at  Princeton  is 
described.  The  general  theme  of  this  effort  is  to  make  the  design  of  VLSI  circuits  as  simi¬ 
lar  to  programming  as  possible.  We  are  trying  to  build  tools  that  do  for  the  VLSI  circuit 
designer  what  the  best  software  tools  do  for  the  implementor  of  large  software  systems 
The  work  described  here  is  a  procedural  language  to  specify  circuit  layouts. 


1.  Introduction 

In  this  paper  we  describe  a  very  important  component  of  any  VLSI 
design  environment:  a  tool  to  automate  the  layout  of  circuits.  This  work 
is  part  of  an  effort  to  create  an  integrated  environment  for  VLSI  design 
(including  layout  systems,  device  and  switch  level  simulators  and  testing 
facilities)  currently  under  way  at  Princeton. 

Our  main  thesis  is  that  the  VLSI  design  task  can  be  profitably  thought 
of  as  a  programming  task,  as  opposed  to  a  geometric  editing  task.  We 
believe  that  much  is  to  be  gained  by  consciously  attempting  to  apply  our 
knowledge  about  programming  to  this  new  activity.  We  have  thus  tried  to 
create  tools  for  the  VLSI  designer  that  incorporate  the  most  useful 
features  of  the  software  develpoment  tools  that  we  are  familiar  with. 

Although  we  feel  that  we  have  had  moderate  success  in  this  endeavor 
we  are  well  aware  of  how  much  room  for  improvement  we  have  left,  and 
would  like  to  help  convince  the  community  of  people  interested  in  the 
design  of  programming  language  and  programming  environments  that 
there  are  fresh  and  important  challenges  in  this  relatively  new  direction. 

A  prototype  of  the  procedural  layout  language  described  in  this 
paper  has  been  operational  for  some  months.  All  figures  given  in  thic 
paper  were  generated  by  the  language  and  all  the  code  fragments  have 
been  used  as  part  of  larger  programs. 

Portion*  of  thi*  paper  appeared  in  the  Proeaodiigf  of  the  IMS  ACM  «ympoo)um  on  Principle*  of  Program 
Ml  language*  and  in  Me  Proceeding*  of  Q»  IMS  DeWgn  AulomaUo  Conference 
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2.  A U:  a  procedural  language  to  describe  layouts 

The  main  feature  of  ALI  as  a  layout  language  is  that  it  allows  its  user 
to  design  layouts  at  a  conceptual  Level  in  which  neither  sizes  nor  posi¬ 
tions  (absolute  or  relative)  of  layout  components  may  be  specified. 
Mostly  as  a  consequence  of  this,  AU  simultaneously  (i)  makes  the  layout 
task  more  like  programming  than  editing,  (ii)  eliminates  the  need  for 
design  rule  checking  after  the  layout  is  generated,  (iii)  permits  the  crea¬ 
tion  of  easy  to  use  cell  libraries  and  (iv)  provides  the  designer  with  the 
mechanisms  to  describe  a  layout  hierarchically  so  that  most  of  the  detail 
at  one  level  of  the  hierarchy  is  truly  hidden  from  all  higher  levels. 

The  notion  of  not  assigning  sizes  or  positions  to  any  object  in  a  layout 
until  the  complete  layout  has  been  described  (similar  to  the  idea  of 
delayed,  binding  in  programming  languages),  sets  ALI  apart  not  only  from 
just  about  all  of  the  graphics  based  layout  editors  we  know  of  ([3],  [4], 
[7],  [13],  [17])  but  also  —  with  the  exception  of  [14]  —  from  most  of  the 
procedural  languages  for  the  layout  task  currently  in  use  or  recently  pro¬ 
posed,  whether  or  not  they  include  a  graphics  interface  ([1],  [4],  [5],  [6], 
[9],  [10],  [15]). 

The  issues  that  we  tried  to  address  with  ALI  are  the  following. 

•  The  creation  of  an  open  ended  tool.  Graphics  editors  tend  to  be  closed 
tools  in  that  it  is  hard  to  automate  the  layout  process  beyond  what  the 
original  design  of  the  system  allowed.  Procedural  languages  are  gen¬ 
erally  much  better  in  this  respect.  However,  the  fact  that  most  such 
languages  require  the  specification  of  absolute  sizes  and  positions, 
makes  the  creation  of  a  general  purpose  library  of  cells  a  hard  task, 
since  information  about  the  sizes  and  positions  of  the  cell  elements  that 
can  interact  with  the  outside  world  has  to  be  apparent  to  the  user  of 
the  library.  The  absence  of  absolute  sizes  and  positions  makes  this 
problem  much  less  severe  in  ALI  The  extensibility  of  ALI  derives  from 
the  fact  that  it  has  been  built  on  top  of  Pascal,  thereby  making  the  full 
power  of  Pascal  available  to  the  designer.  The  generation  of  tools  to 
automate  the  layout  process,  such  as  simple  routers  or  PLA  generators, 
involves  writing  Pascal  routines  to  solve  some  abstract  version  of  the 
problem  and  having  done  so  invoke  ALI  cells  to  generate  the  layouts 

•  Creating  tools  that  are  simple  to  use  and  easy  to  learn.  In  particular, 
we  want  to  avoid  tools  whose  behavior  is  unpredictable.  Many  programs 
which  rely  heavily  on  sophisticated  heuristics  respond  to  small  changes 
in  their  input  with  wholesale  changes  in  their  output.  We  have  main¬ 
tained  a  simple  correspondence  between  the  text  of  an  AU  program  and 
the  resulting  layout  so  that  changes  in  the  layout  can  be  easily  related 
to  changes  in  the  program.  This  decision  has  simplified  the  system  at 
the  cost  of  making  it  less  knowledgeable  about  MOS  circuits 

•  Facilitating  the  division  of  labor.  Large  layouts  have  to  be  produced  by 
more  than  one  designer.  If  the  piece  produced  by  each  designer  is 
specified  in  absolute  positions,  serious  problems  are  likely  to  arise  when 
the  different  pieces  are  put  together,  unless  very  tight  interaction  - 
with  its  attendant  penalties  in  productivity  -  is  maintained  throughout 
the  design.  AU  allows  the  partitioning  of  tasks  in  such  a  way  that  the 
designer  of  a  piece  of  the  layout  does  not  need  to  know  anything  about 
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the  sizes  of  other  pieces  of  the  complete  layout.  For  instance,  on  the 
top  of  fig.  1  three  simple  cells  are  shown  with  the  intended  connections 
between  them  shown  by  dotted  lines;  on  the  bottom  of  the  figure,  the 
pieces  have  been  brought  together  to  form  a  larger  piece.  The  stretch¬ 
ing  that  has  taken  place  has  occurred  without  the  designer  having  to 
plan  for  it  explicitly  while  considering  each  individual  cell. 


RC-l 

Three  separate  cells  and  the  result  of 


connecting  them  along  the  dotted  lines 


i  •  Facilitating  hierarchical  design.  Even  when  a  single  designer  is 

involved,  the  ability  to  view  a  layout  as  a  hierarchy,  with  as  much  infor¬ 
mation  about  lower  levels  completely  hidden  from  higher  levels,  is 
extremely  useful.  In  ALI,  the  information  about  a  given  level  of  the 
hierarchy  needed  at  the  level  immediately  above  is  reduced  by  the 
absence  of  absolute  sizes  and  positions,  to  topological  relations  among 
the  layout  elements  of  the  lower  level  visible  to  the  higher  one. 


•  Reducing  the  life  cycle  cost  of  layouts.  Modifying  a  layout  to  be  fabri¬ 
cated  on  a  new  process,  or  to  make  it  conform  to  a  new  set  of  design 
rules,  is  currently  a  costly  operation.  Yet  successful  designs  seem  to  be 
more  or  less  continuously  updated  as  improved  processes  become  avail¬ 
able  during  their  lifetime.  Fig.  2  shows  two  instances  of  a  simple  layout 
produced  with  ALI.  The  instances  are  the  result  of  running  an  ALI  pro¬ 
gram  twice  changing  exactly  four  constants  in  the  program  in  between 
runs  (those  that  specified  the  sizes  of  power  and  ground  buses).  This 
type  of  flexibility  addresses  the  problem  directly.  An  ALI  program  can 
be  written  naturally  so  that  all  layouts  produced  by  it  are  completely 
free  of  design  rule  violations,  no  matter  what  the  values  of  the  con¬ 
stants  used  in  the  programs.  Therefore  the  need  for  costly  design  rule 
checking  of  different  instances  of  a  layout  (see  fig.  2)  can  be  avoided 
The  same  ALI  program  can  also  generate  layouts  using  different  design 
rules  by  running  it  with  a  new  module  incorporating  the  new  design 
rules. 

•  To  avoid  the  need  /or  special  purpose  computing  equipment.  ALI  can 
be  used  effectively  from  a  standard  ASCII  terminal  in  combination  with 
a  small  plotter  shared  by  several  designers.  All  the  algorithms  used  in 
the  inner  cycle  of  ALI  require  linear  time,  therefore  permitting  the  use 
of  just  about  any  machine  and  guaranteeing  fast  turnaround  on  small 
layouts.  Furthermore  ALI  replaces  design  rule  checking  by  a  hierarchi¬ 
cal  process  that  can  be  performed  separately  on  the  individual  pieces  of 
the  layout.  For  example,  after  checking  that  each  of  the  pieces  shown 
on  the  top  of  fig.  1  is  free  of  design  rule  violations,  their  combination 
shown  on  the  bottom  of  the  same  figure  will  be  guaranteed  by  ALI  to  be 
free  of  rule  violations  regardless  of  the  stretch  that  takes  place  as  a 
consequence  of  connecting  them.  ALI  in  fact  requires  far  fewer  comput¬ 
ing  resources  than  many  design  rule  checking  programs. 

We  feel  that  ALI  succeeds  in  partially  solving  most  of  these  problem 
We  do  not  claim  however  to  have  made  the  layout  task  trivial.  To  use  a 
software  metaphor,  we  feel  that  ALI  elevates  the  work  of  the  layout 
designer  from  absolute  machine  language  programming,  to  programming 
in  a  relocatable  assembler  with  subroutines.  This  not  only  makes  the 
task  more  pleasant  but  makes  new  and  more  powerful  tools  possible 
such  as  loaders,  linkers  and  compilers  in  the  case  of  software.  Similar 
tools  for  the  VLSI  world  -  which  would  indeed  simplify  the  layout  task 
enormously  —  remain,  however,  to  be  written.  ALI  should  stand  or  fall 
with  its  ability  to  allow  such  tools  to  be  built:  whether  we  are  right  in 
believing  that  we  have  a  framework  in  which  these  tools  can  be  more 
easily  implemented  will  not  be  known  until  our  efforts  in  that  direction 
succeed  or  fail. 

The  remainder  of  this  section  is  devoted  to  a  survey  of  the  main 
features  of  ALI  and  a  brief  discussion  of  its  current  implementation 
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Two  ALI  layouts  generated  by  programs 
differing  only  in  the  values  of  four  constants 
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2.1.  An  overview  of  ALI 

The  basic  principles  of  ALI  are  quite  simple.  A  layout  is  regarded  as 
a  collection  of  rectangular  objects  (with  their  sides  oriented  in  the  direc¬ 
tion  of  the  axes  of  a  cartesian  coordinate  system)  and  a  set  of  relations 
among  these  rectangles.  The  ALI  user  specifies  a  layout  by  declaring  the 
rectangles  (also  called  boxes )  of  which  it  is  composed,  and  stating  the 
relations  that  hold  between  them.  ALI  then  generates  a  minimum  area 
layout  that  satisfies  all  the  relations  between  boxes  specified  in  the  pro¬ 
gram.  For  example,  fig.  3  shows  a  trivial  ALI  program  and  the  layout  it 
produces. 

chip  simple 

const 

hwumbtr*  jft- 
length.  =  20; 
width  =  6; 

baxtype 

htype  s  array  [1.. hnumber]  of  metal; 

ear 

i :  integer; 

box 

horizontal :  htype; 
vertical :  metal; 

begin 

fori  :=  1  to  hnumber- 1  do  begin 

above  (  /u>rixonlal[i],  horizontal  [t+  J]  ); 
glueright  (  horizontal[i],  vertical ); 
xmore  (  horizontal  (ij,  length  ) 
end. 

glueright  (  hori*ontal[hmim6rr] ,  vertical ); 
xmore  ( Aorixontal[/inumber],  length  ); 
xmore  (  vertical,  width  ) 

and. 


ng.3 

A  simple  ALI  program  and  the  layout  it  produces 


This  program  looks  very  much  like  a  Pascal  program:  it  consists  of  a 
declarative  part,  followed  by  an  executable  part.  To  declare  a  box  the 
user  specifies  its  name  ( horizontal  or  vertical  in  the  example),  and  if  s 
type,  ( metal  -  a  predefined  type  -  in  the  example).  The  standard  box 
types  correspond  to  the  layers  of  the  physical  layout.  As  the  example 
also  shows,  the  ALI  user  can  define  structured  objects  (an  array  in  the 
example).  Further  details  on  the  type  structure  of  ALI  can  be  found  in 
section  2.2.1. 

The  relations  between  the  rectangles  that  make  up  a  layout  are 
specified  in  ALI  through  calls  to  a  small  set  of  primitive  operations  in  the 
executable  part.  All  such  operations  take  as  arguments  boxes  and  possi¬ 
bly  values  of  standard  Pascal  types  (integers  in  our  example).  In  our 
example  above,  glueright  and  xmore  are  primitive  operations.  The  primi¬ 
tive  above  specifies  that  its  first  argument  must  appear  above  the  second 
one  in  the  final  layout,  the  primitive  glueright  extends  its  first  argument 
to  the  right  \'i  intersect  ts  second  argument,  and  xmore  makes  the  size 
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of  its  first  arguments  along  the  x  axis  at  least  as  large  as  the  value  of  the 
second  argument.  Note  that  in  this  example  ALI  has  determined  the 
minimum  separation  between  the  horizontal  elements  as  well  as  the 
minimum  sizes  of  boxes  not  specified  by  zmore  (such  as  the  height  of  the 
horizontal  metal  lines)  by  accessing  a  table  of  design  rules.  More  infor¬ 
mation  about  the  type  structure  and  the  primitive  operations  of  ALI  is 
given  in  the  next  section. 

When  an  ALI  program  is  executed  it  generates  two  kinds  of  informa¬ 
tion.  It  produces  a  set  of  linear  inequalities  involving  the  coordinates  of 
the  corners  of  the  boxes  in  the  layout  as  variables.  These  inequalities, 
which  embody  the  relations  between  the  rectangles  of  the  layout,  are 
then  solved  to  generate  the  positions  and  sizes  of  the  layout  elements  A 
brief  description  of  the  problems  involved  in  this  step  can  be  found  in 
section  2.3.2.  The  program  also  produces  connectivity  information  about 
the  rectangles  in  the  layout.  11118  information  can  then  be  used  by  a 
switch  level  simulator  that  predicts  the  behavior  of  the  circuit  as  laid  out 
whithout  having  to  perform  the  ususal  "node  extraction"  analysis. 

In  order  for  the  layouts  produced  by  an  ALI  program  to  be  free  of 
design  rules,  the  program  must  be  complete,  in  that  every  pair  of  rectan¬ 
gles  in  it  must  be  related  in  some  way.  Two  rectangles  may  be  related 
explicitly  in  the  user  program  by  virtue  of  being  arguments  to  a  primitive 
operation,  or  they  may  be  related  through  the  transitivity  of  the  separa¬ 
tions.  The  reason  for  this  strong  requirement  is  to  prevent  the  area 
minimization  process  from  shoving  together  rectangles  that  were 
intended  to  be  separate  (see  section  2.3.3  fo  a  discussion  of  complete¬ 
ness). 

ALI  helps  the  designer  to  achieve  this  goal  by  generating  certain  rela¬ 
tions  between  layout  elements  in  an  automatic  fashion,  and  by  checking 
on  request  whether  this  condition  is  satisfied.  It  is  however  the  responsi¬ 
bility  of  the  user  to  make  an  ALI  program  complete  in  this  sense,  as  the 
computational  cost  of  doing  any  sophisticated  inference  (beyond  the 
transitivity  of  relations  such  as  above)  is  prohibitive  [16]. 


2.2.  Main  features  of  ALI 

This  section  describes  how  ALI  appears  to  its  user  Its  three  subsec¬ 
tions  deal,  in  turn,  with  the  type  structure,  the  primitive  operations  of 
the  language  and  the  cell  mechanism.  ALI  has  been  built  on  top  cf  Pascal 
and  has  inherited  most  of  its  features.  In  the  interest  of  shortening  this 
section  we  have  assumed  a  certain  familiarity  with  the  general  features  of 
Pascal. 


2.2.1.  Type  structure 

As  the  example  of  fig.  3  shows,  the  objects  manipulated  by  ALI  are 
declared  by  stating  their  name  and  their  type.  The  types  of  ALI  have  the 
same  structure  as  the  Pascal  types.  Objects  can  be  of  a  simple  type 
(boxes)  or  of  a  structured  type. 
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There  are  a  small  number  of  standard  types,  all  of  them  simple.  The 
standard  types  correspond  to  the  layers  of  the  process  to  be  used  to 
fabricate  the  layout  (metal,  poly,  diff,  impl,  cut  and  glass  in  the  NMOS 
version  currently  implemented)  plus  the  type  virtual,  used  to  name 
bounding  boxes  and  having  no  physical  reality  in  the  fabricated  circuit. 
For  example,  in  the  program  of  fig.  2,  the  declaration 

vertical  :  metal 

specifies  that  the  rectangle  named  vertical  on  the  final  layout  should  be 
on  the  metal  layer.  ALI  will  use  this  information  to  generate  constraints 
on  its  minimum  size  and  its  separation  from  other  layout  elements. 

Structured  types  are  of  two  flavors:  array  (a  collection  of  objects  of 
the  same  type)  and  bits  (a  collection  of  objects  of  heterogeneous  types, 
much  like  records  in  Pascal),  which  correspond  directly  to  the  array  and 
record  structured  types  of  Pascal.  ALI,  like  Pascal,  permits  the  creation 
of  new  user  defined  types  that  can  be  either  simple  or  structured  For 
example,  in  fig.  3,  the  fragment 

htype  =  array  [l  h.num.ber]  of  metal 

inside  the  boxtype  section  of  the  program,  creates  a  new  type,  htype, 
each  object  of  that  type  made  up  of  a  number  of  metal  rectangles,  and 
the  fragment 

horizontal  :  htype 

inside  the  box  section,  creates  an  object  of  that  type  named  horizontal 

In  a  similar  fashion  the  type  declaration 

shiftbus  =  bus 

phi,  ph2 :  metal; 
vdd  :  metal; 
data  :  diff, 
gnd  :  metal 

end 

creates  a  user  defined  type,  allowing  the  user  to  create  objects  which 
consist  of  four  metal  boxes  and  a  diffusion  box.  The  types  of  the  com¬ 
ponents  of  structured  types  are  arbitrary:  the  user  can  define  arrays  of 
buses,  or  buses  containing  arrays. 

The  accessing  of  the  elements  of  arrays  and  buses  is  done  as  in  Pas¬ 
cal.  Thus  if  *  is  of  type  htype  then  x[i]  refers  to  the  i-th  element  of  x.  and 
if  y  is  of  type  shiftbus  then  y.data  refers  to  the  diffusion  box  of  y 

Although  the  structured  objects  are  generally  used  by  ALI  simply  as  a 
naming  mechanism,  they  are  also  used  in  conjunction  with  the  ceil 
mechanism  (discussed  in  section  2.2.3)  to  automatically  generate  separa¬ 
tions  between  boxes.  We  will  be  more  precise  on  this  point  when  we 
describe  the  cell  mechanism  of  ALI. 
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Uke  Pascal,  AL1  is  a  strongly  typed  language.  The  primitive  opera¬ 
tions  know  about  certain  type  restrictions  and  generate  type  mismatch 
errors  if  operations  are  attempted  with  rectangles  of  inappropriate 
types. 


2.2.2.  Primitive  operations 

The  relations  between  the  rectangles  that  make  up  a  layout  are 
specified  in  ALI  through  calls  to  a  small  set  of  primitive  operations.  All 
such  operations  take  boxes  (i.e.,  objects  of  simple  types)  as  arguments 
In  the  program  of  fig.  3,  above,  glueright  and  imore  are  primitive  opera¬ 
tions. 

It  is  not  important  to  know  the  actual  primitive  operations  of  the 
current  version  of  ALI  to  understand  its  operation.  As  a  gross  measure  of 
its  complexity  we  can  say  that  the  system  currently  implemented  - 
based  on  NMOS  as  described  in  [12]  —  has  about  twenty  primitive  opera¬ 
tions  which  can  be  arranged  in  the  following  groups: 

1  Separation  primitives:  such  as  above  in  fig.  3.  which  specify  that  their 
arguments  must  be  separated  in  a  certain  direction  in  the  final  layout. 
The  minimum  amount  of  space  between  boxes  separated  in  this  manner 
depends  on  their  types  and  is  supplied  by  ALI  from  a  table  of  design 
rules. 

2  Connection  primitives:  such  as  glueright  in  fig.  3,  to  specify  that  their 
arguments  --  which  must  be  boxes  in  the  same  layer  --  are  to  be  joined 
in  a  particular  manner. 

3  An  inclusion  primitive,  inside,  that  specifies  that  one  box  is  to  be 
placed  inside  another.  The  minimum  distances  between  their  edges  are 
again  suplied  by  ALI  from  a  table  of  design  rules. 

4  Minimum  size  primitives:  such  as  zmore  in  fig.  3,  which  specify  the 
minimum  size  of  a  box  along  a  certain  direction.  Default  minimum 
sizes  are  provided  by  ALI  from  a  design  rule  table. 

5  Transistor  primitives,  which  create  depletion  mode  and  pass  transis¬ 
tors 

6  Contact  Primitives,  which  create  contacts  between  layers  and  connect 
boxes  to  them 

Note  that  no  absolute  positions  or  dimensions  for  any  rectangle  can 
be  specified  with  these  primitives.  All  the  rectangles  of  a  layout  can  be 
stretched  and  compressed  (up  to  a  minimum  size)  and  all  can  float  in  any 
direction.  If  one  single  characteristic  is  to  be  used  to  separate  ALI  from 
other  layout  systems,  this  must  be  it.  Most  of  the  power  of  Al.i  and  most 
of  the  problems  one  faces  in  its  implementation  are  consequences  of  this 
fact. 

It  is  important  to  remember  that  in  order  for  a  layout  produced  by 
ALI  to  be  free  of  design  rule  violations,  any  two  rectangles  in  it  must  be 
related  in  some  way.  ALI  will  make  no  inferences  as  to  the  relations 
between  boxes  beyond  those  implied  by  the  transitivity  of  some  primitive 
operations  (i.e.,  if  above  (a,  b)  and  above  (b,  c)  are  stated,  above  (a  c) 
need  not  be  stated).  Although  the  system  generates  a  good  number  of 


relations  automatically  for  the  user,  particularly  in  connection  with  the 
cell  mechanism  (see  the  next  subsection),  there  is  still  a  fair  amount  of 
drudgery  left  for  the  user  in  making  sure  that  this  requirement  is  met  A 
brief  discussion  on  the  computational  complexity  of  the  automatic  gen¬ 
eration  of  relations  between  boxes  can  be  found  in  section  2.3.3. 


2.2.3.  Cells 

Perhaps  the  most  powerful  feature  of  ALI  is  its  procedure-like 
mechanism  for  the  definition  and  creation  of  cells.  A  cell  is  a  collection 
of  related  rectangles  enclosed  in  a  rectangular  area.  Rectangles  that  are 
inside  a  cell  are  of  two  types:  Local  which  are  invisible  to  the  outside,  or 
parameters  which  can  interact  in  a  simple  and  well  defined  manner  with 
rectangles  outside  the  cell. 

A  cell  is  defined  by  specifying  its  local  objects,  its  formal  parameters 
and  the  relations  among  all  of  them.  Once  a  cell  has  been  defined,  it  can 
be  instantiated  as  many  times  as  desired  by  specifying  the  actual  param¬ 
eters  for  the  instance,  much  the  same  way  as  one  invokes  a  procedure  or 
function  in  a  procedural  language.  The  result  of  instantiating  a  cell  is  to 
create  a  brand  new  copy  of  the  prototype  described  in  the  cell  definition 
with  the  formal  parameters  connected  to  the  actual  parameters. 

A  cell  definition  is  made  up  of  a  header,  in  which  the  formal  parame¬ 
ters  are  described,  a  set  of  local  box  declarations  and  a  body  in  which  the 
relationship  between  the  parameters  and  the  local  boxes,  as  well  as  those 
among  local  boxes,  are  specified. 

The  header  describes  the  names  and  types  of  the  parameters  and  the 
side  of  the  bounding  rectangle  through  which  they  come  into  contact  with 
the  inside  of  the  cell.  The  header  of  a  cell  (using  the  type  shiftbus 
defined  in  section  2.2. 1)  and  an  instance  of  it  are  shown  in  fig  4 


call  shift  (Mil:  shifttms;  right  r  :  shiftbus  ) 


t  W  -it  -'*-**♦'<*-'•'**  B 


ng.4 

A  sample  cell  definition  header 
and  an  instance  of  the  cell  defined 
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Cells  may  have  any  number  of  parameters  on  each  of  their  four  sides. 
The  order  in  which  they  are  listed  in  the  cell  header  describes  their  rela¬ 
tive  positions.  Horizontal  parameters  (i.e. ,  those  touching  the  cell  on  the 
left  or  right)  are  assumed  to  be  listed  in  top  to  bottom  order  and  vertical 
parameters  in  left  to  right  order. 

The  body  of  a  cell  is  very  much  like  an  ALI  program.  For  example,  flg. 
5  shows  a  complete  cell  definition  that  consists  of  a  variable  number  of 
shift  cell  instances  connected  sequentially  together  with  two  of  its 
instances.  Note  that  cells  are  instantiated  by  the  create  statement,  and 
that  the  parameter  list  of  the  cell  contains  both  box  parameters  and 
other  parameters  (an  integer  in  this  case)  in  separate  lists.  Note  also 
that  recursion  has  been  used  to  define  this  cell;  this  highlights  the  fact 
that  ALI  has  the  full  power  of  Pascal  at  its  disposal. 

When  an  instance  of  a  cell  is  created  it  can  be  given  a  name,  provided 
that  the  name  given  has  been  declared  as  a  rectangle  of  the  standard 
simple  type  virtual.  The  relationship  of  the  rectangle  bounding  a  newly 
created  cell  to  any  other  rectangle  of  the  layout  can  be  specified  in  the 
standard  manner  by  calls  to  the  primitive  operations.  This  is  a  vital 
feature  since  in  many  cases  (i.e.,  above,  below...)  stating  a  relation 
between  two  cell  instances  c,  and  ce  immediately  implies  a  relation 
between  every  pair  of  rectangles  r,  and  rg  such  that  r,  is  part  of  c,  and  r2 
part  of  c8. 

There  are  two  important  ways  in  which  the  cell  mechanism  helps  in 
the  automatic  generation  of  constraints  between  boxes.  When  an  object 
of  a  structured  type  is  passed  as  a  parameter  to  a  cell,  its  component 
boxes  are  separated  from  top  to  bottom  (if  it  is  a  left  or  right  argument) 
or  from  left  to  right  (if  it  is  a  top  or  bottom  argument).  The  order  of  the 
separation  is  determined  by  applying  recursively  the  following  rules: 
array  elements  are  separated  ordered  by  their  indices  and  bus  elements 
in  the  order  in  which  they  were  specified  in  the  bus  declaration.  Thus,  in 
the  example  of  fig.  5,  the  components  of  parameter  inbus  would  be 
separated  from  top  to  bottom.  The  second  mechanism  involves  the 
automatic  separation  of  cells  that  share  a  parameter;  thus  in  the  exam¬ 
ple  of  fig  5,  the  individual  instances  of  shift  are  separated  automatically, 
since  adjacent  instances  share  a  parameter. 

The  cell  mechanism  gives  the  ALI  user  the  ability  to  describe  layouts 
in  a  truly  hierarchical  manner.  A  proper  ALI  design,  very  much  like  a  well 
structured  program,  will  consist  of  a  hierarchy  of  cell  instances  with  only 
a  small  amount  of  information  at  a  given  level  (the  parameters  of  the  cell 
instances  at  that  level)  being  visible  from  the  immediately  higher  level. 
For  example,  the  layout  given  in  fig.  2  consists  of  four  instances  of  the 
same  cell  stacked  vertically.  That  cell  in  turn  is  defined  in  terms  of  three 
other  cells,  one  of  them  being  the  cell  shown  in  fig.  1,  which  is  in  turn 
defined  in  terms  of  three  other  cells. 

Much  of  the  power  and  generality  of  the  cell  mechanism  of  ALI  comes 
from  the  absence  of  absolute  positions  and  sizes  in  a  layout  specification. 
In  particular,  two  instances  of  the  same  cell  may  have  radically  different 
sizes  depending  on  the  actual  parameters  used  to  create  them,  as 
exemplified  by  figs.  1,  2  and  5.  We  believe  that  no  cell  mechanism  can  be 
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cell  shiftregister  (  left  tabus  ;  shift  bus; 
right  outbus  :  shxftbus  ) 

( length  :  integer  ); 

box 

temp  :  shift  bus; 

begin 

if  length  =  1  then 

create  shift  ( tabus,  out  bus  ) 

dee  begin 

create  shift  ( tabus,  temp  ); 
create  shiftregister  ( temp,  outbus  )  ( length  -1 ) 
and 
anti 


ilc.5 

A  cell  definition  and  two  instances  of  it 
generated  by  a  simple  ALI  program 


said  to  be  truly  general  unless  the  sizes  of  its  arguments  and  local  rec¬ 
tangles,  as  well  as  their  relative  distances  are  determined  at  the  tune  the 
cell  is  instantiated. 

There  are  some  penalties  involved  in  the  use  of  the  cell  mechanism 
In  particular,  ALI  generates  separations  between  cells  in  a  mariner  which 
is  oblivious  to  what  is  inside  them.  That  is,  the  minimum  separation 
between  cells  as  far  as  ALI  is  concerned,  is  the  maximum  of  all  the 
minimum  separations  for  two  layers  in  the  design  rules,  thus  creating  a 
certain  wastage.  We  believe  that  this  penalty  will  be  generally  a  small 
percentage  of  the  total  area  and  is  well  worth  the  advantages  gained  by 
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the  ability  to  separate  cell  instances  as  units. 

Another  source  of  wastage  is  the  fact  that  cells  are  restricted  to  be 
bounded  by  a  rectangle,  so  the  packing  of  cells  that  have  irregular  shapes 
results  in  a  certain  amount  of  unused  space.  The  rectangular  shape  of 
the  cells  is  a  fundamental  characteristic  of  ALI:  The  introduction  of  irreg¬ 
ularly  shaped  cells  is  simply  not  possible  without  completely  redesigning 
the  language  and.  However,  the  waste  introduced  because  of  this  restric¬ 
tion  can  be  avoided  in  most  particular  cases  through  some  code 
modifications. 


2.3.  Implementation  issues 

The  previous  section  described  the  user  view  of  ALL  In  this  section 
we  discuss  briefly  some  of  the  problems  to  be  solved  when  trying  to  go 
from  an  ALI  program  to  a  layout  that  satisfies  the  relations  stated  in  it 
We  first  give  an  overall  description  of  the  system  as  currently  imple¬ 
mented,  then  discuss  the  method  used  to  assign  locations  and  sizes  to 
the  layout  elements  and  then  the  concept  of  completeness  and  how  it  is 
checked. 


2.3.1.  Overall  implementation 

The  current  version  of  our  system  has  been  implemented  as  follows. 
The  ALI  program  is  first  translated  into  standard  Pascal.  The  resulting 
Pascal  program  is  then  compiled  and  linked  with  a  precompiled  set  of 
procedures  that  implement  the  primitive  operations  and  the  resulting 
object  module  is  then  run.  The  output  of  this  program  (generated 
entirely  by  the  primitive  operations)  is  a  set  of  linear  inequalities  and 
connectivity  relations  among  the  layout  elements.  The  inequalities  are 
then  solved  to  generate  a  layout  or  examined  by  a  program  that  checks 
their  logical  completeness,  and  the  connectivity  information  can  be  used 
to  simulate  the  circuit  laid  out. 

The  design  rules  are  incorporated  as  a  table  which  is  used  by  the 
primitive  operations  to  produce  the  linear  inequalities.  Thus  changing 
the  design  rifles  for  our  system  requires  only  to  change  this  table 


2.3.2.  Placement 

As  explained  above,  one  of  the  results  of  running  an  ALI  program  is  a 
set  of  linear  inequalities  that  embody  the  relations  between  the  layout 
elements.  These  inequalities  are  of  the  following  simple  form: 

xt-xj3td  (dfcO) 

where  the  variables  are  the  coordinates  of  the  corners  of  the  boxes  that 
form  the  layout  and  the  constants  are  either  user  supplied  (as  in  the 
second  argument  of  the  xmore  primitive,  for  instance)  or  extracted  from 
the  table  of  design  rules  by  the  system  itself. 
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This  set  of  inequalities  should  be  solved  so  as  to  generate  placements 
for  the  boxes  that  compose  the  layout  in  such  a  way  as  to  minimize  iLs 
total  area.  In  order  to  perform  this  task  efficiently,  we  require  that  no 
inequality  in  the  set  involve  both  x  and  y  coordinates.  This  restriction 
allows  us  to  minimize  the  total  area  by  minimizing  the  maximum  x  and  y 
coordinates  of  any  point  independently,  at  the  cost  of  reducing  the  range 
of  the  relations  between  boxes  that  we  can  express.  We  cannot,  for 
instance,  handle  rectangles  whose  sides  are  not  parallel  to  the  cartesian 
axes  or  express  aspect  ratios  directly. 

We  have  now  a  sufficiently  simple  problem  so  it  can  be  solved  in  time 
proportional  to  the  number  of  inequalities  in  our  set.  (All  layouts  that 
can  be  expressed  in  AU  can  be  generated  by  a  program  that  produces  a 
constant  number  of  inequalities  per  rectangle).  This  is  done  by  a  version 
of  the  topological  sort  algorithm  [11]  applied  to  the  x  and  y  coordinates 
independently.  This  algorithm  assigns  to  each  point  the  lowest  possible 
coordinate  while  minimizing  the  largest  coordinate  of  all  points. 

The  form  of  the  inequalities  that  we  allow  is  rather  restrictive;  it  is 
sufficient  however,  to  describe  the  design  rules  given  in  [12]  for  NMOS, 
and  the  efficiency  gained  in  return  for  this  simplicity  seem  to  us  like  a 
good  tradeoff.  A  more  subtle  consequence  of  the  simplicity  of  the  ine¬ 
qualities  and  the  method  we  use  to  solve  them  is  that  undesirable 
stretching  can  occur,  since  we  have  no  way  to  specify  a  maximum  size  for 
any  object.  This  is  not  a  common  occurrence  and  the  user  can  in  all 
cases  guard  against  such  stretching  by  the  careful  selection  of  the  primi¬ 
tive  operations  used.  It  is  nonetheless  an  additional  burden  placed  on  the 
designer. 

The  choice  of  an  efficient  placement  algorithm  over  expressibility 
power  and  a  reduced  degree  of  user  convenience  has  been  quite  cons¬ 
cious  in  this  particular  case.  We  feel  that  every  reasonable  measure 
should  be  taken  to  keep  the  complexity  of  the  placement  problem  linear, 
given  that  the  size  of  layouts  is  currently  large  (1QT  rectangles)  and  is 
growing  fast.  Widening  the  class  of  linear  inequalities  acceptable  is 
almost  certain  to  make  linear  time  solutions  impossible  [2]. 


2.3.3.  Completeness 

AU  programs  do  not  involve  absolute  sizes  or  positions  of  boxes,  and 
are,  to  a  great  extent, independent  of  the  design  rules.  These  characteris¬ 
tics  make  it  clearly  desirable  to  insure  that  the  layout  described  by  a 
program  will  be  free  of  design  rule  violations  in  a  way  other  than  check¬ 
ing  the  finished  layout.  The  following  paragraphs  describe  a  way  of  insur¬ 
ing  freedom  from  design  rule  violations  in  a  manner  that  is  independent 
of  the  actual  design  rules  used  to  generate  the  Anal  placement.  The 
description  may  be  somewhat  cryptic;  the  interested  reader  is  referred 
to  [16]  for  further  details. 

A  layout  generated  by  an  AU  program  is  complete  if  for  any  two 
boxes  a  and  b  whose  types  make  it  possible  for  them  to  interact  in  the 
final  layout,  either 
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(i)  a  and  b  are  explicitly  stated  to  be  in  contact  by  some  primitive 
operation,  or 

(ii)  a  and  6  are,  explicitly  or  through  the  transitivity  of  primitive  rela¬ 
tions,  stated  to  be  separated  in  either  the  *  or  the  y  direction  by  a 
minimum  amount  which  depends  on  their  types. 

From  this  definition,  it  should  be  clear  that  testing  completeness  of  a 
cell  instance  involves  computing  the  transitive  closure  of  a  graph.  There¬ 
fore  the  complexity  of  the  operation  will  be  0(ns),  where  n  is  the  number 
of  boxes  in  the  cell.  It  is  thus  not  feasible  to  test  a  large  layout  for  com¬ 
pleteness  in  a  direct  way. 

Fortunately,  completeness  can  be  checked  hierarchically.  Let  us 
look  only  at  the  objects  at  the  highest  level  of  the  hierarchy  of  boxes  that 
defines  a  layout,  i.e.,  those  boxes  (including  cell  boundaries)  defined  glo¬ 
bally  in  the  ALI  program  that  generated  the  layout.  If  these  objects  are 
related  in  a  complete  manner  and  the  cell  instances  used  at  this  level  are 
also  complete,  then  the  whole  layout  is  complete. 

Thus  one  can  check  the  completeness  of  a  layout  by  successively 
checking  cell  instances  for  completeness,  thereby  reducing  the  complex¬ 
ity  of  the  process  to  0(m3)  where  m  is  the  largest  number  of  boxes  local 
to  a  cell  instance  in  the  layout.  This  process  can  be  reduced  further, 
since  not  every  cell  instance  needs  to  be  checked.  For  example,  if  a  cell 
is  defined  by  a  straight  line  program,  checking  one  instance  for  com¬ 
pleteness  suffices,  as  one  instance  of  the  cell  will  be  complete  if  and  only 
if  all  of  its  instances  are  [16].  The  case  of  cells  with  branches  and  itera¬ 
tion  is  not  quite  as  simple.  Yet  we  are  confident  --  and  our  experience 
tends  to  confirm  this  belief  --  that  checking  the  completeness  of  a  few 
carefully  selected  instances  of  any  cell  definition  will  be  enough  to 
guarantee  that  the  cell  definition  is  complete. 

The  end  result  is  that  completeness  has  the  flavor  of  a  static,  almost 
syntactic,  property  for  all  non  malicious  examples,  and  is  much  easier  to 
check  in  a  well  structured  layout  than  design  rule  freedom  by  the  stan¬ 
dard  means  on  the  final  layout. 

Finally,  a  word  about  the  possibility  of  taking  an  incomplete  layout 
specification  and  automatically  completing  it.  The  general  problem  of 
generating  an  optimal  completion  is  NP-Complete,  but  the  simpler  ver¬ 
sion  of  generating  any  completion  for  graphs  embedded  in  a  grid  (as  our 
layouts  are)  seems  to  be  solvable  in  0(n*)  steps.  The  question  of  how 
much  area  will  be  wasted  by  such  a  completion  algorithm  will  have  to  wait 
for  some  experimentation,  but  there  is  no  question  of  its  usefulness. 


2.4.  Experience  with  ALI 

The  current  implementation  of  ALI  has  shown  the  soundness  of  most 
of  our  original  ideas.  The  system  is  efficient  and  the  language  easy  to 
learn,  and  the  layout  it  produces  are  relatively  dense  (for  example,  an  ALI 
program  written  without  concern  for  area  optimization  produced  a  layout 
which  was  about  thirty  per  cent  larger  than  a  similar  layout  packed  by 
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hand  on  a  graphics  editor.  Unfortunately,  this  evidence  has  been  gath¬ 
ered  mostly  from  people  who  had  a  hand  in  designing  or  implementing 
ALI.  Perhaps  a  more  reliable  evaluation  of  ALI,  ought  to  wait  until  a  sub¬ 
stantial  number  of  users  not  involved  in  its  design  can  give  an  informed 
opinion.  We  hope  to  obtain  this  evidence  before  long,  since  ALI  is 
currently  being  used  in  a  VLSI  design  course. 

The  fact  that  very  little  effort  was  invested  in  error  recovery  for  the 
sake  of  expediency  in  getting  a  prototype  running,  and  that  no  mechan¬ 
ism  for  integrating  separately  produced  layout  pieces  was  provided  make 
the  current  system  useful  mostly  for  teaching  purposes  and  experimen¬ 
tation.  It  must  be  emphasized  that  this  is  a  result  of  implementation 
choices,  and  not  of  any  intrinsic  limitation  on  the  approach  we  have 
taken. 

The  problems  of  the  current  system  which  we  plan  to  address  with 
the  next  version  are  the  following: 

(1)  Memory  requirements.  The  solution  of  the  system  of  linear  inequali¬ 
ties  requires  large  amounts  of  memory.  We  will  use  a  different  algo¬ 
rithm  which  is  slightly  less  efficient  in  terms  of  time  but  requires  an 
order  of  magnitude  fewer  memory  locations  for  a  typical  large  lay¬ 
out. 

(2)  Pascal  problems.  The  current  ALI  has  exactly  the  same  type  struc¬ 
ture  as  Pascal.  The  lack  of  generic  types  and  dynamic  arrays  has 
made  the  task  of  writing  general  purpose  tools  (PLA  generators, 
routers...)  inside  ALI  more  difficult  than  it  ought  to  be.  The  next  ALI 
will  have  the  notions  of  generic  types  and  dynamic  arrays 

(3)  Connecting  primitives  Certain  objects,  such  as  contacts,  are  used 
frequently  enough  to  warrant  making  them  part  of  the  language. 

(4)  Separate  “compilation"  facilities.  Clearly,  large  layouts  will  have  to 
be  generated  in  pieces,  which  is  something  that  our  current  system 
cannot  do. 
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1.  Introduction 

The  problem  we  address  in  this  paper  is  an  embedding  problem  for  a  class 
of  graphs  which  we  call  rectilinear  graphs.  These  graphs  are  important  in  many 
VLSI  layout  problems.  In  fact,  this  problem  arose  in  the  implementation  of  ALI 
[6,7],  a  procedural  language  for  VLSI  design  currently  under  development  at 
Princeton.  An  embedding  algorithm  can  be  used  to  automate  the  production  of 
VLSI  layouts  in  many  procedural  design  systems. 

Consider  the  following  model  for  VLSI  Layout  design.  A  VLSI  Layout  is 
described  hierarchically  using  cells  and  wires  that  connect  the  cells  together. 
Each  cell  C  is  enclosed  within  a  rectangle  /?(C),  and  has  four  tuples  of  pins,  one 
each  for  the  left,  top.  right,  and  bottom  of  rectangle  R{C).  Each  wire  w  is 
denoted  by  a  pair  of  pins  (PiPj),  such  that  ft  and  Pj  are  pins  of  different  cells 
and  are  of  opposite  types.  For  example,  if  ft  is  a  right  pin  then  Pj  should  be  a 
left  pin.  Given  such  a  description  of  a  VLSI  layout,  our  aim  Is  to  produce  an 
embedding  of  the  description  on  the  plane,  such  that  (i)  no  two  bounding  rec¬ 
tangles  touch  each  other,  (ii)  the  pins  appear  in  the  correct  order  on  the 
bounding  rectangles.  (OL)  the  wires  are  straight  and  rectilinear,  and  (t*/)  no  two 
wires  cross  each  other.  Later  on,  we  can  fill  each  bounding  rectangle  J 7(C)  with 
the  embedding  of  the  cell  C  in  the  same  manner. 

The  restriction  that  wires  cannot  be  bent  may  seem  unrealistic,  but  this  is 
certainly  the  case  in  many  design  systems  including  ALI.  If  a  wire  has  to  be 
bent,  the  user  specifies  that  by  breaking  up  the  wire  into  several  straight  vires 
and  placing  cells  at  each  of  the  turn  points  of  the  wire.  In  ALI,  for  example,  the 
user  can  incorporate  routing  algorithms  in  a  ALI  program  to  determine  how  the 
wires  are  to  be  bent.  The  restriction  that  wires  cannot  cross  implies  that  we  are 
dealing  with  the  wires  on  a  single  layer.  For  a  layout  with  multiple  layers,  it  is 
clearly  necessary  that  the  wires  on  each  layer  do  not  oross. 

To  solve  the  above  embedding  problem,  it  is  enough  to  consider  only  a  sim¬ 
ple  restriction,  where  each  bounding  rectangle  has  at  most  one  pin  on  each  side. 
Ve  can  then  treat  the  bounding  rectangles  as  vertices  and  the  wires  as  edges, 
which  leave  in  one  of  the  four  cardinal  directions.  We  give  the  name  rectilinear 
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graphs  to  such  graphs.  The  embedding  problem  of  rectilinear  graphs  is  our 
main  concern  in  this  paper. 

For  VLSI  applications,  we  need  efficient  algorithms  to  recognize  and  then 
actually  embed  rectilinear  graphs.  In  this  paper,  we  present  an  0(n)  recogni¬ 
tion  algorithm  and  an  0(n*)  embedding  algorithm,  where  n  is  the  number  of 
vertices  in  the  graph.  Thus,  a  hierarchically  described  VLSI  layout  with  cell 

instances  Ci,C8,  •  •  •  ,C^  can  embedded  in  time  «(£%■>.  where  n,  is  the 

<*I 

number  of  pins  in  cell  instance  Q . 

An  embedding  of  a  rectilinear  graph  is  just  a  relative  placement  of  the  ver¬ 
tices  (  cells  )  on  a  rectangular  grid,  such  that  no  two  edges  cross.  Some  of  the 
relative  placement  information  is  already  present  in  the  description  of  a  rectil¬ 
inear  graph.  For  example,  if  (a.b)  is  the  rightgoing  edge  of  vertex  a,  then  a 
should  be  to  the  right  of  b,  and  a.b  should  be  on  the  same  horizontal  grid  line. 
Hence,  an  embedding  can  be  viewed  as  a  "completion"  of  the  rectilinear  graph 
description.  We  showed  in  a  different  paper  [8]  that  the  completion  problem  for 
a  slightly  more  relaxed  VLSI  layout  model  is  NP-complete.  In  light  of  this  result, 
the  results  in  this  paper  have  become  more  important. 

In  section  2.  we  present  formal  definitions  of  rectilinear  graphs  and  their 
embeddings.  In  section  3.  we  mention  some  properties  of  rectilinear  graphs.  We 
discuss  some  topological  properties  of  the  embeddings  in  section  4.  A  necessary 
and  sufficient  condition  for  biconnected  rectilinear  graphs  to  be  embeddable  is 
presented  in  section  5.  A  similar  condition  for  arbitrary  rectilinear  graphs  is  the 
main  result  In  section  6.  We  also  describe  a  0(n)  recognition  algorithm  in  this 
section.  In  section  7,  we  use  the  ideas  of  the  previous  sections  to  obtain  an 
0(nz)  embedding  algorithm  An  Important  subclass  of  rectilinear  graphs  is  dis¬ 
cussed  in  section  B.  In  section  8,  we  discuss  extensions  and  open  problems.  For 
definitions  of  graph  theoretic  terminology  used  in  this  paper,  please  refer  to 
[1.2]- 

2.  Definition  of  the  Problem 

First  we  give  a  formal  definition  of  a  rectilinear  graph. 

Definition  2.1:  A  rectilinear  graph  G  is  a  triple  [V.  E,  A),  where  V  is  the  vertex 
set,  E  is  the  edge  set,  and 

A  :VxV-»  E  \j  }ej  ,  where  E  =  | L,R,D,U\ 

is  a  vertex  ordering  relation  with  the  following  properties 
for  every  a  ,b  ,e  €  V and  X  €  E 

(t)  A((a,b))  =  e  «♦  (a,b  j  £  E 

(ordering  is  specified  only  between  adjacent  vertices) 

(«)  A((a,b))  =  L  •  A((b,a))  =  R,  A((a,b ))  =  D  «•  A((b .a))  s  U. 
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(iii)  X((a,b))  =  X  •*  X((c,b))  *  X.  Vc  ?  a  (no  overlapping  edges). 

Each  vertex  in  a  rectilinear  graph  has  degree  at  most  four,  and  each  edge 
(a.b),  as  it  goes  from  one  vertex  a  to  the  another  b ,  has  a  nonempty  label  on  it. 
which  in  the  embedding  will  indicate  the  direction  (left,right,  down,  or  up)  in 
which  the  edge  leaves  vertex  a.  There  can  be  at  most  one  edge  with  a  particular 
label  emanating  from  each  vertex.  The  undirected  graph  G(V,E)  will  be 
referred  to  as  the  underlying  graph..  Figure  2. 1  (like  all  other  figures)  gives  an 
illustration  of  a  rectilinear  graph. 

— *  R 


Figure  2  1 
A  rectilinear  graph 


Now  we  define  what  sort  of  an  embedding  we  are  looking  for. 

Definition  2.2:  An  embedding  of  a  rectilinear  graph  G(K.f’.X)  on  a  rectangular 
grid  is  given  by  two  mappings  x,  y :  V-*  Z  (  the  integers  )  which  are  the  x  and  y 
coordinates  respectively  of  the  vertices.  These  mappings  obey. 

1.  the  ordering  relation,  X,  i.e.  for  all  edges  (a.b  J  e  E 
X((a.b))  =  L  -*  y(a)=y(6),  *(a)>x(b), 

X((a,b))  =  /?  -•  y(a)  =  y(b),  x(a)<x(b), 

X((e.b))  =  D  -*  x(a)  =  x(b),  y(a)>y(b), 

X((o,b))=£7  -*  x(a)=x(6),  y(a)<y(b). 

2  Planarity,  no  two  edges  cross,  i.e.  for  each  pair  of  non-adjacent  edges 
(a.b  j.  (c  ,dj  such  that  X((a,b))  =  R  and  \((c , d ))  =  U,  the  condition 

*(a)  <  x(c)  <  x(b)  and  y(c)sy(a)sy(i) 

does  not  hold. 

An  embedding  of  a  rectilinear  graph  on  a  rectangular  grid  is  one  in  which 
the  vertices  are  placed  at  grid  points,  the  edges  run  along  grid  lines  in  the  direc¬ 
tions  given  by  their  labels,  and  no  two  edges  cross  each  other  except  if  they 
share  a  vertex  We  say  that  a  rectilinear  graph  is  embeddable  if  it  has  an 
embedding  We  will  show  in  the  next  section  that  not  all  rectilinear  graphs  are 
embeddable. 
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Now  our  main  problem  can  be  stated  simply:  Given  a  rectilinear  graph 
G(  V,E,X).  is  it  embeddable,  and  if  yes,  find  an  embedding. 

3.  Some  Coomenti  on  Rectilinear  Graphs 


• - — •  *  i 

Figure  3.1 

Tiro  nonembeddable  rectilinear  cycles 


Figure  3.2 

Two  nonembeddable  rectilinear  graphs 
whose  biconnected  components  are 

In  this  section  we  list  some  properties  of  rectilinear  graphs  and  their 
embeddings  These  will  give  an  indication  of  why  the  problem  is  not  trivial  and 
why  it  is  different  from  other  embedding  problems,  and  in  particular,  planar 
graph  embedding  [3,5], 

1.  Embeddabilxty  is  a  hereditary  property.  Subgraphs  are  defined  in  the  usual 
fashion,  but  here  the  labels  of  edges  are  inherited  This  is  obvious,  but 
worth  mentioning,  because  this  will  be  used  in  the  proofs. 

2  If  each  connected  component  of  a  rectilinear  graph  is  embeddable  then  the 
graph  itself  is  embeddable.  So,  without  loss  of  generality  we  will  restrict 
ourselves  to  connected  rectilinear  graphs 

3  Rectilinear  graphs  with  nonplanar  underlying  graphs  are  clearly  not 
embeddable.  So  it  is  not  interesting  to  consider  those  graphs  However,  not 
every  rectilinear  graph  with  a  planar  underlying  graph  is  embeddable.  In 
figure  3  1 ,  we  have  two  simple  cycles  which  are  not  embeddable. 
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4.  In  contrast  with  planarity,  embeddability  is  not  a  property  determined  by 
the  biconnected  components.  Figure  3.2  provides  an  illustration  of  this 
fact. 

5  This  problem  is  a  restriction  of  an  NP-complete  problem  [0,11].  For  each 
wire  ui.  we  are  given  its  orientation  (horizontal  or  vertical),  and  a  set  Vw  of 
vertices.  Hie  wire  w  has  to  pass  through  each  vertex  in  the  set  Vw  (the  ver¬ 
tices  could  be  touched  in  any  order).  Then,  the  embedding  problem 
becomes  NP-complete. 

6.  If  we  relax  the  rectilinearity  of  the  edges  and  impose  only  the  cyclic  order¬ 
ing  of  the  edges  at  each  vertex,  then  there  is  an  0(  |  V|)  algorithm  [10].  The 
cyclic  orderings  automatically  determines  the  faces  of  the  embedding  (if 
one  exists).  Thus  a  embeddable  rectilinear  graph  has  a  unique  embedding  in 
this  sense. 


b 

\({abcdefgh))  =  LURDRVL 
Figure  4  1 

The  extension  of  A  to  paths 

There  is  a  natural  way  to  extend  the  function  X  to  paths  and  cycles  in  the 
graph  as  follows.  Given  a  path  P  =  (v0,vi,  ••  •  ,vt)  we  define 
A(P)  =  X((t/0.v1))X((i'i,vg))  •  •  ■  X((vj_i,W|))  We  define  a  similar  extension  for 
cycles  where  now  vt  =  v0.  A  becomes  a  mapping  that  associate*  '■'tth  each  path 
or  cycle  in  the  graph  a  string  in  E*  which  is  the  concatenation  of  labels  along  the 
path  or  cycle.  Note  that  strings  containing  RL,  DU,  LR,  UD  as  substrings  do 
not  represent  paths.  Also  the  direction  in  which  we  traverse  a  path  and  the 
starting  point  in  a  cycle  are  important.  An  example  of  this  mapping  can  be  found 
in  figure  4.3. 

Next  we  define  two  topological  actions  on  rectilinear  graphs  These  actions 
will  simplify  a  rectilinear  graph  while  preserving  Its  topological  structure.  Let  G 
be  a  rectilinear  graph 

Action  1  -  Edge  Contraction:  Let  (abed)  be  a  path  in  G  such  that  both  6  and  c 
have  degree  2,  and  A((a6cd))  =  XYX  where  X.Y  e  E.  Contract  the  edge  (6,c)  to 
the  vertex  b.  The  resulting  path  (abd)  will  have  A((abd))  =  XX  We  abbreviate 


I 


-6- 


Figure  4.2 

Edge  contraction  and  vertex  deletion 
this  action  by  XYX  -*  XX  (  figure  4.2(1) ). 

Action  2  -  Vertex  Deletion:  Let  (obc )  be  a  path  in  G  such  that  vertex  b  has 
degree  2,  and  X((obc))  =  XX  where  X  £  £.  Delete  the  vertex  b  and  introduce 
the  edge  (a,c).  The  resulting  edge  (a,c)  will  have  X((a ,c))  =  X.  We  abbreviate 
this  action  by  XX  -*  X  (  figure  4.2(2) ). 

In  a  natural  way  we  can  define  inverses  for  the  above  two  actions  which  we 
will  refer  to  as  edge  expansion  and  vertex  addition  respectively. 

Lemma  4.1:  Let  G  be  a  rectilinear  graph  and  G'  be  the  graph  resulting  from  G 
by  the  application  of  a  sequence  of  the  above  four  actions  Then  G'  is  also  rectil¬ 
inear  and  moreover  G'  is  embeddable  if  and  only  if  G  is  embeddable. 

Proof:  The  proof  is  easy  and  is  left  to  the  reader.  • 


RDLDR  -  R 


RDLDR  ULDLURDL V  -  RDLV 


Figure  4.3 

Simplification  of  a  path  and  a  cycle 


Definition  4.1:  Given  a  string  7CE*  representing  a  path  or  a  cycle,  the 
timplifUtd  form  ’y  of  7  is  obtained  by  repeatedly  applying  the  reduction  rules 
XYX  ■*  XX  and  XX  -*  X,  where  X,Y  e  E,  until  they  cannot  be  applied  any  more. 
If  7  represents  a  cycle  then  it  is  treated  as  a  cyclic  string. 

In  figure  4.3  we  give  a  path  and  a  cycle  along  with  their  simplified  forms. 
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Lemma  4.2:  Every  string  y  £  2*  has  a  unique  simplified  form 

Proof:  The  replacement  system  defined  by  the  two  reduction  rules  have  the 

Church-Rosser  property  [B].  • 

Definition  4.2:  A  square  is  one  of  the  cyclic  strings  LURD  or  LDRV 

Sometimes  we  may  distinguish  between  two  squares  by  their  starting  labels 

Definition  4.3:  A  spiral  is  a  path  which  cannot  be  simplified  Equivalently  a 
spiral  is  a  substring  of  (LURD)*  or  ( LDRU )*. 

Lemma  4.3:  Every  path  is  embeddable. 

Proof:  Every  spiral  is  embeddable.  Since  any  path  simplifies  to  a  spiral  (by 
definition ).  by  lemma  4.1  it  is  also  embeddable.  • 

So  it  is  the  cycles  which  make  the  problem  nontrivial  The  following  lemma 
is  a  crucial  fact  about  cycles 

Lemma  4.4:  A  cycle  is  embeddable  if  and  only  if  it  simplifies  to  a  square. 

Proof:  if:  A  square  is  embeddable  and  hence  by  lemma  4.1  any  cycle  which 
simplifies  to  a  square  is  also  embeddable. 

only  if:  Let  /  be  an  embeddable  cycle  and  \(/)  =  y.  By  lemma  4  1,  the 
cycle  defined  by  y  is  also  embeddable.  Let  jy|  =  n.  Look  at  the  embedding  of  y 
Since  it  has  no  crossings  the  embedding  is  a  simple  polygon.  Therefore  the  inte¬ 
rior  angle  of  this  polygon  sum  to  (n-2).l80*.  Since  y  is  a  spiral  all  its  interior 
angles  are  90°.  The  only  solution  to  n.  90  =  (n-2)lBD  is  n  =  4.  Therefore  ^  is  a 
square.  • 

The  proof  of  the  previous  lemma  suggests  another  useful  characterization 
of  embeddable  cycles.  Going  along  a  cycle  /  =  vtvg  •  •  •  vnvt  in  the  counter¬ 
clockwise  direction,  let  us  denote  by  p/(v*)  the  angle  at  vertex  vx,  which  is  the 
angle  between  (i u-iWi)  and  (vj,v1+1). 


Lemma  4.5:  A  cycle  /  =  Vji/g  •  •  •  n  se  4  is  embeddable  if  and  only  if 

*(/)  =  £*/("<)  =  (n± 2)180'  . 

t*i 

Proof:  Suppose  f  is  embeddable,  then  its  embedding  is  a  simple  polygon. 
Depending  on  whether  we  sum  the  Interior  angles  or  exterior  angles  we  should 
get  (n  ±2).1B0°. 

To  prove  the  sufficient  part  we  show  by  induction  of  n  that  /  simplifies  to  a 
square.  The  possible  values  for  ?/(«<)  are  90,180,270°.  The  basis  for  induction  is 
n  =  4.  In  this  case  the  given  sum  of  the  angles  is  either  360°  or  1060°,  which 
implies  that  each  angle  is  either  90°  or  360°  respectively.  So  /  must  be  a 
square  by  itself. 

Assume  that  the  claim  is  true  for  all  values  less  than  n  and  let  n  >  4  If  for 
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some  i,  f/(u<)  =  260°.  then  A(t/<_|i/(v*«.i)  =  XX  We  can  apply  vertex  deletion  at 
v4  to  obtain  /'  *  vjv*  •  •  •  Vi-iVt+i  ■  •  •  v»v,.  Then 

fit')  =  fit)  ~  ff(vt)  =  ((n-l)±2).l80°.  and  by  induction  we  are  done. 

This  leaves  the  case  where  all  angles  are  either  90°  or  270°.  Since  n  >  4 
and  p(/  )  =  (n±2).180°  not  all  the  angles  can  be  equal.  Hence  there  must  be  a  * 
such  that  f/(vg)  ff/(vk4.x).  Hence  we  have  X(t>fe_,vkvk4.tvfe+S)  =  XYX  Apply 
edge  contraction  to  obtain  /'  =  vt  •  •  •  v*_iut/k+z  ■  ■  •  vnvy  The  edge  contrac¬ 
tion  removed  360°  from  the  angle  sum  and  added  180°.  Hence 
fit  )  -  ((*i-l)±2).l80“.  • 

Definition  4.4:  A  complement  of  a  path  P  with  respect  to  a  square  a  is  any 
path  Pcin  the  graph  such  that  PPe is  a  cycle  which  simplifies  to  \{PPe)  =  a. 

Lemma  4.6:  Given  a  path  P,  all  its  complements  with  respect  to  a  square  a. 
which  have  the  same  start  and  end  labels,  have  a  unique  simplified  form. 

Proof:  Let  A (P)  =  a  =  XxXg  2*  Since  a  is  a  spiral  we  have  X+  -  Xx  for 
i  =  j  (4).  Assume  that  *  >  4  and  that  the  spiral  a  and  the  square  a  are  either 
both  clockwise  or  both  counterclockwise.  Then  a  must  be  a  substring  of  a.  Since 
o  is  a  cyclic  string  we  can  assume  that  a  =  XxXtXsX 4. 

Let  Pc  be  a  complement  of  P  with  respect  to  a  and  let  A(Pt)  =  0.  Since 
k  >  4.  0  must  spiral  in  the  opposite  direction  to  a.  Since  both  a  and  0  are 
simplified,  a0  can  be  simplified  only  at  the  borders  between  the  two  strings. 
Write  0  =  01020a.  such  that  03a0i  =  a  We  are  allowed  to  shift  03  because  cl0  is  a 
cyclic  string.  Then  it  is  clear  that  0X  e  (Ak.9A'k,JI«,e{  and  09  e  02  is 

the  'essential  part'  of  0.  Since  |a|  =  Jfc  and  0  -  4,  we  must  have  \0t  |  =  *-4. 
From  the  possible  values  of  0X  and  03.  and  the  fact  that  0  is  a  spiral  opposite  in 
direction  to  a  ,  we  can  conclude  that  0e  =  XA  We  used  k  >  4  in 

order  for  0Z  not  to  be  an  empty  string.  Therefore 
0  =  [t,Xic.Xk.3Xk]Xk^ ,  •  •  •  X4UJCx.XxX4],  which  is  unique  but  for  the  start  and 
end  labels.  The  arguments  in  the  cases  where  a  and  a  are  in  opposite  directions 
and  for  k  s  4  are  similar.  • 

S.  Bi connected  Rectilinear  Graphs 

In  this  section  we  discuss  an  algorithm  for  recognizing  biconnected  rectil¬ 
inear  graphs.  Note  that  the  ordering  relation  A  induces  a  cyclic  ordering  of  the 
edges  incident  at  each  vertex  v.  For  convenience  we  will  need  the  following 
definition. 

Definition  5.1:  Let  v  be  a  vertex  in  a  rectilinear  graph  G.  Define  Lc(v)  to  be  the 
cyclic  list  of  the  neighbors  of  v  in  G  in  the  counterclockwise  order. 

Using  these  lists,  we  can  define  the  essential  notion  of  a  candidate  face  of  a 
biconnected  rectilinear  graph. 
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C/’1(t)  =  2, 1,5,4, 3, 2  and  C/*(0  =  1.2, 7, 6.5.4, 3, 2 
Figure  51 
Candidate  faces 

Definition  5.2:  Let  G  =  ( V.E.k)  be  a  biconnected  rectilinear  graph.  With  each 
edge  e  =  (t'l.vg),  Ju,  >  we  associate  two  lists  of  vertices  called  candidate 
facts  CF  i(e)  and  CFs(e )  which  are  defined  as  follows. 

CF ,(0)  =  1/,,1/g.  •  •  ■  ,wk.w*+i  where  i/t  ¥■  for  ij  *  w*+i,  i  ¥■  j  and  =  vt 

for  some  i,  1  <i  <  Jfc-1,  such  that  for  each  l,  1  <  l  <4  +  1,  is  the  successor 
of  V|_,  in  the  cychc  list  Lc(vt).  CFt(s)  is  similarly  defined  but  starting  with 

Wg.O,. 

It  is  easy  to  see  that  CFt  and  CFt  are  uniquely  defined  An  illustration  of 
this  definition  is  given  in  figure  5.1. 

We  now  need  a  lemma  about  biconnected  undirected  graphs.  Let  us  define 
a  biconnected  graph  to  be  minimal  if  for  every  edge  #  in  the  graph  G-*  is  not 
biconnected.  The  following  lemma  is  taken  from  [2]  and  is  stated  without  proof. 

Lemma  5.1:  If  G  is  a  minimal  biconnected  graph  having  at  least  four  vertices 
then  G  contains  a  vertex  of  degree  two. 

Lemma  5.2:  In  any  biconnected  graph  G  which  is  not  a  simple  cycle,  there  is  a 
simple  path  P  =  (vl,vg),[v2,va),  ■  •  •  .(1/,,-i.tv),  r  a  2,  with  the  intermediate  ver¬ 
tices  (  If  any  )  vit  i  *  l,r  all  having  degree  2,  such  that  the  graph  G'  =  G-P  is 
biconnected. 

Proof:  Transform  the  given  graph  G  to  another  graph  G"  by  replacing  all  paths 
of  the  form  P  -  (vi,v2),(ve,t;s),  •  •  •  .(t/y-j.w,.)  where  the  vertices  «<.  i  *  l.r  all 
have  degree  2,  by  the  edge  {vx,vT).  So  for  each  edge  •  in  G"  we  have  a 
corresponding  path  P,  in  G.  Note  that  the  degree  of  any  vertex  in  G"  is  at  least 
three.  If  G"  has  multiple  edges  between  some  two  vertices,  say  u  and  su,  then  In 
G  there  must  be  at  least  two  parallel  paths  between  u  and  w .  Since  G  is  not  a 
simple  cycle  any  one  of  those  paths  will  serve  our  purpose.  If  G"  does  not  have 
multiple  edges  then  it  must  have  at  least  4  vertices.  By  lemma  5.1  G"  cannot  be 
minimal  Therefore  there  is  an  edge  e  in  G"  such  that  G"- s  is  biconnected. 

t  For  convenience  we  assume  that  V  is  a  aet  of  integer*. 
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which  implies  that  G-P,  is  also  biconnected.  ■ 

4  The  following  theorem  gives  a  necessary  and  sufficient  condition  for  a  bicon¬ 

nected  rectilinear  graph  to  be  embeddable. 


M.i  \lf  tv* 

Figure  5.3 

The  two  cycles  / ,,/  e  and  the  path  P 


Theorem  5.1:  Let  G  =  ( K,£\A)  be  an  biconnected  rectilinear  graph  with  at  least 
three  edges.  Then  G  is  embeddable  if  and  only  if  for  each  edge  e  in  the  graph 
both  the  candidate  faces  CF,(e)  and  CFe(a )  represent  simple  embeddable 
cycles  in  the  graph  (i.e.  the  starting  and  ending  vertices  are  identical,  and  it 
simplifies  to  a  square).  Moreover,  If  the  graph  is  embeddable  each  such  distinct 
I  candidate  face  corresponds  to  a  face  in  the  planar  embedding. 

Proof: 

only  if:  Supposing  for  some  edge  e  =(vj,vg),  CFj(e)  is  not  a  cycle,  i.e 
CF i(e)  =  Vi,vs,  ■  •  ■  1  with  =  t/*+i  for  some  i,  1  <  <  <  *-1.  Suppose  G 

is  embeddable.  Look  at  the  cycle  w*.  •  •  •  in  the  embedding  Suppose 

that  the  edge  («t_|, vt)  is  inside  this  cycle.  There  can  be  no  other  edges 
{u.vj),  is  i  inride  this  cycle,  otherwise  u  would  have  appeared  instead  of 
vj+i  in  CFi(e)  From  the  planarity  of  the  embedding,  there  An  be  no  path  from 
v(_i  to  Vi  other  than  the  edge  (v4_i,Vi).  This  contradicts  the  biconnectedness  of 
G.  The  case  where  (v(_i,vt)  is  outside  the  cycle  is  similar  (  both  cases  are  dep- 
'  icted  in  figure  5.2 ). 


Suppose  CF j(«)  is  a  cycle  but  is  not  embeddable.  Since  CFx(t )  is  a  sub¬ 
graph  of  G,  G  itself  cannot  be  embeddable.  Similar  arguments  holds  for  GFg(e). 

if.  The  proof  of  this  part  is  by  induction  on  the  number  of  edges.  The  basis 
for  the  induction  are  simple  embeddable  cycles,  for  which  the  claim  is  true  by 
lemma  4.4.  Assume  that  the  claim  is  true  for  any  biconnected  rectilinear  graph 
which  has  less  than  k  edges.  Let  G  be  a  biconnected  rectilinear  graph  which  is 
not  a  simple  cycle  and  which  has  k  edges.  By  lemma  5.2,  there  is  a  simple  path 
P  =  (v1.va),(v8,vs)<  •  •  •  .(ty-i.ty)  with  the  vertices  vt,  i  w  l,r  all  having  degree 
2.  such  that  the  graph  G'  =  G-P  is  biconnected.  i»,  and  iv  will  have  degree 
greater  than  two.  Also  assume  that  >  v*  and  e,s  =  (v  ,  ,v  £) . 

Since  all  our  candidate  faces  are  cycles,  if  an  edge  a  lies  on  a  candidate 
face  /  then  either  CF i(e)  =  /  or  CF e(s )  =  / .  So  each  edge  will  be  present  in 
exactly  two  of  these  candidate  faces.  Hence  the  path  P  will  appear  in  CF x(*ig) 
and  its  reverse  path  will  appear  in  CFs(e  12).  Let 
/ 1  =  CFx{k12)  =  vx.v2.  •  •  •  .vr.u1(  •  •  •  .Uj.Vi  . 

/ g  =  CFa(*u)  =  vr.vT-x,  •  •  •  .«,.«>!.  •  •  •  ,wk.Vr  .  and 

/ 8  =  >rMX.  ■  ,«*, VX. 

It  follows  from  the  definition  of  the  candidate  faces  /  4  md  /z  that  the  vertices 
tiy.i/j.ui,  appear  consecutively  in  that  order  in  Lc{vx)  and  that  uv.iv_i.it, 
appear  similarly  in  Le(vr)  (  see  figure  5.3  ).  Therefore  for  each  edge  in  / ,  or  /  z 
which  is  not  in  P,  the  new  candidate  face  in  G'  will  be  /  3  which  is  a  simple  cycle. 

We  still  have  to  show  that  / a  is  embeddable.  Since  /,  and  /t  are  both 
embeddable  p(/ ,)  =  (r  +y  ±2).  180°  and  p(/g)  =  (r+k  ±2).  180°  .  However,  since 
/,  and  fa  share  the  edge  e,s  it  is  implied  by  the  definition  of  candidate  faces 
that  «*(/ 1)  =  (r+^+  2). 160°  and  f{fg)  -  (r+*+2).160°  is  impossible.  With  a  little 
bit  of  algebraic  manipulation  we  can  show  that  f[f  8)  =  {{j  +A+2)±2).180*.  Since 
f3  has  j  +  k +2  vertices  by  lemma  4.5,  it  is  embeddable  Thus  the  candidate  faces 
for  G'  are  the  same  as  those  for  G,  excepting  for  /9  replacing  the  two  faces  / , 
and  /  g.  So  for  each  edge  in  G‘  its  two  candidate  faces  are  again  simple  embedd¬ 
able  cycles.  By  induction  hypothesis  G‘  is  embeddable  and  each  distinct  candi¬ 
date  face  corresponds  to  a  face  in  its  embedding  The  orderings  of  the  edges  at 
the  vertices  v,  and  iv  imply  that  the  end  edges  (v,.vg)  and  (iv-,.iv)  of  the  path 
P  are  both  trying  to  go  inside  the  face  corresponding  to  /s. 

We  are  left  to  show  that  we  can  add  the  path  P  back  without  destroying 
embeddability.  Find  any  rectilinear  path  P'  in  the  face  corresponding  to  /s  in 
the  embedding  of  G'.  that  starts  and  ends  with  A((v,,vg))  and  A((vr_,.vr)) 
respectively.  This  is  clearly  possible  although  we  may  have  to  extend  the  grid  in 
order  for  P'  to  lie  on  the  grid  lines.  P'  creates  a  face  in  the  embedding  with  the 
path  Pi-  «vU|Ug  -  *  -  Ujt/,.  If  /  s  is  not  the  outside  face  then 
\{PXP)  -  Kfi)  -  A(/*,  P)  -  o.  The  case  when  /  3  is  the  outside  face  is  slightly 
more  complicated.  There  are  two  such  different  paths  P'  depending  on  the  new 
outer  face  that  is  created.  However,  for  one  of  the  two  the  above  holds  and  sup¬ 
pose  this  is  the  one  we  chose.  By  definition  both  P  and  P'  are  complements  of 
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P\  with  reaped  to  a.  they  also  share  the  same  start  and  end  labels,  and  by 
lemma  4.6  we  have  \[P)  -  Therefore  G'+P‘  can  be  obtained  from  G  by 

applying  a  sequence  of  the  four  topological  actions,  and  since  G'+P'  is  embedd¬ 
able,  by  lemma  4.1  G  is  also  embeddable.  It  is  easy  to  see  that  the  two  new 
faces  we  get  after  inserting  P  in  the  embedding  of  G  correspond  to  / }  and  /  8 .  • 
The  above  theorem  leads  to  the  following  algorithm  for  recognizing  embedd¬ 
able  biconnected  rectilinear  graphs.  The  algorithm  also  outputs  the  faces  of  the 
embedding  if  the  graph  happens  to  be  embeddable 

Algorithm  chack-biconnscted(G), 
begin 

If  G  is  an  edge  then  return; 
if  \E\  >  3 j  K |  — 6  then 
begin 

write  (  Hot  embeddable'), 
quit 

end 

for  each  edge  e  do 
begin 

mark[e  ,l]:=  false, 
mark[e  .2].=  false 

end, 

for  each  edge  e  do 
for  i;=  1  to  2  do 
begin 

if  not  mark' e  ,t]  then 
begin 

/  :=  candidate-face  f  e,  i  ), 

If  not  embed  -cycle  (/ )  then 
begin 

write (  Hot  embeddable' ), 
quit 

end; 

for  each  edge  e'  =  (v,,v2)  in  f  do 
if  Vj  >  vz  then  mark[s  ',l]:=  true 
else  mark[e  ',2].s=  true; 
output  (f  ) 

end 

end 

end. 

Boolean  function  tmbed-cycle(f )  returns  value  true  if  /  is  an  embeddable 
cycle.  If  /  is  a  cycle  then  we  simplify  using  the  reduction  rules  and  check  if  we 
end  up  with  a  square.  This  can  be  done  in  time  linear  in  the  size  of  / .  Function 
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call  candidate-/ ace (e ,i)  returns  the  candidate  face  CF<(e)  and  the  function  can 
be  implemented  exactly  as  described  in  definition  5.2.  In  the  calls  to  this  func¬ 
tion,  each  edge  e  can  be  traversed  at  most  twice,  due  to  the  flags  mark[e  ,1]  and 
mark[e  ,2],  Therefore  the  algorithm  runs  in  time  0{  j  V  |).  We  conclude  this  sec¬ 
tion  with  a  lemma  which  will  let  us  identifiy  the  outer  face  in  a  rectilinear  graph. 


!«■«»  5.3:  Let  G  be  an  embeddable  biconnected  rectilinear  graph.  For  all 
interior  faces  /  in  the  embedding  of  G,  9(f)  -  (n-2).l80“,  and  for  the  unique 
exterior  face  /, ,  ?(/,)  =  (n+2).180*. 

Plroof:  Consider  the  embedding  of  G.  The  faces  of  the  embedding  are  deter¬ 
mined  by  G.  and  are  simple  polygons  in  the  plane.  By  the  definition  of  9.  for 
every  interior  we  count  the  interior  angles,  and  for  the  exterior  face  we  count 
the  exterior  angles  The  lemma  follows.  (  Remember  that  if  G  is  a  simple  cycle, 
the  embedding  has  two  faces  ).  ■ 
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Figure  5.4 
Shapes  U  and  W 


Lemma  5.4:  Let  G  be  an  embeddable  biconnected  rectilinear  graph,  /,  the 
exterior  face  in  its  embedding  and  v  a  vertex  on  /, .  If  9/,(v)  =  180°  then  G 
can  be  embedded  inside  a  polygon  of  shape  U,  as  shown  in  figure  5.4a.  If 
9/m  =  270"  then  G  can  be  embedded  inside  a  polygon  of  shape  W,  as  shown  in 
figure  5  4b. 

Proof:  Easy  and  left  to  the  reader.  ■ 


6.  Articulation  Vertices 

In  this  section  we  examine  the  conditions  under  which  the  embeddability  of 
the  biconnected  components  of  the  graph  Implies  the  embeddability  of  the 
graph  itself.  Clearly,  this  will  depend  on  the  way  components  meet  at  articula¬ 
tion  vertices.  In  figure  3.2,  we  showed  two  examples  of  nonembeddable  rectil¬ 
inear  graphs,  each  of  which  decomposes  into  two  embeddable  biconnected  rec¬ 
tilinear  graphs. 

In  those  cases,  the  two  biconnected  components  are  not  "compatible"  at 
the  articulation  vertex.  However,  the  situation  need  not  be  so  local.  Figure  6.1 
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Figure  6.1 

Decompositions  of  nonembeddable  graphs 


depicts  two  nonembeddable  graphs,  each  of  which  decomposes  into  three 
embeddable  biconnected  components,  so  that  the  components  meeting  at  each 
articulation  vertex  are  compatible.  Note  that  an  edge  is  a  (trivial)  biconnected 
component. 

If  v  is  an  articulation  vertex  in  a  graph  G,  then  its  removal  results  in 
several  connected  subgraphs  G*  of  G  We  will  refer  to  the  subgraphs  ta+f.  as 
the  subgraphs  meeting  at  v .  Throughout  this  section  we  will  implicitly  assume 
that  we  are  dealing  with  rectilinear  graphs  whose  biconnected  components  are 
embeddable 

Definition  8.1:  Let  £?,  and  B%  be  two  nontrivial  biconnected  components  of  a 
rectilinear  graph  G  that  share  an  articulation  vertex  v.  Then  Bx  and  Bz  are  said 
to  interlace  if  the  horizontal  edges  at  v  belong  to  B\  and  the  vertical  edges 
belong  to  Bz  We  also  say  that  v  is  an  interlace  vertex.  Any  articulation  vertex 
that  does  not  have  this  property  is  said  to  be  intrrlace-free. 

Lemma  6. 1:  A  rectilinear  graph  G  which  has  an  interlace  articulation  vertex  v  is 
not  embeddable. 

I¥oof:  Let  Bv  and  Bz  be  the  two  biconnected  components  sharing  the  vertex  v 
Since  Bx  and  Bz  are  nontrivial,  the  horizontal  edges  at  v  lie  on  a  cycle  in  B j  and 
the  vertical  edges  lie  on  a  cycle  in  Bz  It  is  impossible  to  draw  G  on  the  plane 
without  these  two  cycles  crossing.  ■ 

Definition  6.2:  Let  Bx  and  Bz  be  two  non  interlacing  biconnected  components  of 
G  that  share  an  articulation  vertex  v,  and  assume  B\  is  nontrivial.  Then  Bz  is 
said  to  be  inside  Bx  (  or  Bx  dominates  Bg  at  v  )  if  either  (i)  v  is  not  on  the  exte¬ 
rior  face  of  Bx,  or  (ii)  edges  (v,u)  and  ( v,w )  at  the  vertex  v  are  on  the  exterior 
face  of  Bx,  and  u,i u  are  consecutive  in  that  order  in  Ac(v)  (  note  that  they  are 
always  consecutive  in  Lg^v)  ).  If  neither  Bx  dominates  Bz  nor  Bz  dominates  B j, 
then  B\  and  Bz  are  said  to  be  outside  each  other. 
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The  intuition  behind  the  above  definition  is  that  in  the  embedding,  one 
i  biconnected  component  must  lie  wholly  inside  some  face  of  the  other  if  one 

edge  of  it  does.  This  is  due  to  the  planarity  criterion.  Clearly,  if  biconnected 
components  B  i  and  Bg  that  share  an  articulation  vertex  v  dominate  each  other, 
the  graph  is  not  embeddable  ( this  is  the  case  in  figure  3.2b  ). 

Let  B i  and  Bg  be  two  biconnected  components  of  a  graph  G  that  share  an 
l  articulation  vertex  v ,  such  that  Bx  dominates  B g.  Let  G'  be  the  subgraph  of  G 

meeting  at  v,  that  contains  Bg.  If  G  is  embeddable  then  in  any  embedding  of  G, 
all  of  G'  should  lie  inside  one  face  of  Bx.  This  suggests  extending  the  relation 
"dominate"  as  follows. 

Definition  6.3:  Let  B  =  \Bx,Bg,  •  -  *  ,Bm\  be  the  set  of  biconnected  components 
of  G.  We  say  that  Bx  dominates  Bj  if  there  exists  a  biconnected  component  Bk 
and  an  articulation  vertex  v ,  such  that  (i)  B*  and  Bi  share  v ,  (ii)  Bx  dominates 
Bt  at  v ,  and  (tit)  B}  and  Bk  are  both  subgraphs  of  the  same  connected  sub¬ 
graph  meeting  at  v . 

Let  us  denote  by  V(G)  the  vertex  set  of  the  graph  G  and  by  E{G)  the  edge 
set 

1-gmma  6.2:  If  in  a  rectilinear  graph  G.  there  exists  some  pair  of  biconnected 
components  Bx  and  Bg  that  dominate  each  other,  then  G  is  not  embeddable. 
Proof:  If  Bx  and  Bg  share  an  articulation  vertex  v ,  then  as  mentioned  earlier  G 
is  not  embeddable.  Suppose  that  B \  and  Bg  are  disjoint.  Since  Bx  and  Bg  dom¬ 
inate  each  other,  there  must  be  articulation  vertices  vx,ve.  biconnected  com¬ 
ponents  Bx,Bg,  and  subgraphs  Gx,Gg,  such  that  for  i  =  1,2,  (i)  Bx  and  Bx'  share 
Vx,  (ii)  Bi  dominates  Bx  at  vx,  and  (tit)  Gj  is  one  of  the  subgraphs  meeting  at  v< 
and  contains  Bi'.  Let  us  assume  that  G  is  embeddable.  FYom  (i)  e  V(C,). 
(»)  Gx  lies  wholly  inside  Bx  in  the  embedding,  and  (iii)  V{Gx)r,V{B x)  =  HM.  we 
can  conclude  that  ve  must  be  properly  inside  a  polygon  defined  by  the  face  / 1  of 
Bx  containing  vx.  Similarly  v ,  should  be  properly  inside  the  polygon  defined  by  a 
face  fgotBg  containing  vz.  Therefore  some  vertices  of  f  g  must  lie  outside  / 1 
f  and  the  two  faces  must  intersect,  and  hence  G  is  not  embeddable.  • 

Given  a  rectilinear  graph  G,  with  set  of  biconnected  components  B  and  set 
of  articulation  vertices  A.  we  can  construct  a  tree  T  of  biconnected  components 
such  that  V(T)  =  AuB,  and  E{T)  =  |(v,B)  |  v  e  A  B  e  B.  v  e  V{B)\. 

1  6.3:  Let  G  be  a  rectilinear  graph  with  the  set  of  biconnected  com¬ 

ponents  B  and  tree  of  biconnected  components  T.  Let  B  be  a  leaf  in  the  tree  T 
which  is  adjacent  to  an  articulation  vertex  v  of  degree  2  in  T.  If  B  dominates  B' 
the  other  biconnected  component  adjacent  to  v  in  T,  then  B  dominates  every 
other  biconnected  component  in  B. 

Proof:  The  only  two  subgraphs  meeting  at  v  are  B  and  G-B+v  and  the  proof 
follows  from  definition  6.3.  • 
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If  no  two  biconnected  components  dominate  each  other,  then  the  relation 
"dominate"  induces  a  partial  order  on  B.  A  nondominating  element  in  this  partial 
order  is  a  biconnected  component  which  does  not  dominate  any  biconnected 
component. 

Corollary  6.1:  If  for  a  rectilinear  graph  G,  "dominate"  is  a  partial  order,  then 
there  exists  a  nondominating  biconnected  component  which  is  a  »eaf  in  the  tree 
T  of  biconnected  components 

Proof:  Any  trivial  biconnected  component  (  which  is  just  an  edge  )  must  be  non¬ 
dominating.  If  any  vertex  in  T  (  corresponding  to  an  articulation  vertex  in  G  )  is 
adjacent  to  two  leaves,  then  either  the  two  leaves  are  nontrivial  and  not  dom¬ 
inating.  or  one  of  them  is  a  trivial  biconnected  component.  If  no  vertex  in  T  is 
adjacent  to  two  leaves,  then  all  leaves  are  adjacent  to  vertices  of  degree  2.  and 
there  are  at  least  two  such  leaves.  If  two  of  these  leaves  are  dominating,  then  by 
lemma  6.3  the  two  leaves  dominate  each  other  which  is  a  contradiction  that 
"dominate"  is  a  partial  order.  In  fact  all  of  these  leaves  must  be  nondominating.  ■ 

Theorem  6.1:  Let  G  be  a  rectilinear  graph  and  B  its  set  of  biconnected  com¬ 
ponents.  G  is  embeddable  if  and  only  if 

(i )  every  biconnected  component  B  in  B  is  embeddable. 

(ii)  every  articulation  vertex  in  G  is  interlace-free,  and 
(Hi)  "dominate"  induces  a  partial  order  onB 

Proof:  The  necessary  part  follows  from  lemma  6  1  and  lemma  6.2. 

The  sufficient  part  is  shown  by  induction  on  the  number  of  vertices.  The 
basis  for  induction  is  any  biconnected  rectilinear  graph  Let  G  be  not  bicon¬ 
nected  with  |K(G)|  =  n.  Assume  that  the  claim  is  true  for  all  smaller  graphs. 
Look  at  the  tree  T  of  biconnected  components.  By  corollary  8.1,  there  exists  a 
leaf  B  in  T  which  is  nondominating.  Let  v  be  the  articulation  vertex  shared  by  B 
and  G'  -  G-B+v.  the  rest  of  the  graph.  G'  being  a  subgraph  of  G  also  satisfies 
the  conditions  of  the  claim.  By  induction  hypothesis  G'  is  embeddable.  By  condi¬ 
tion  (i).  B  is  also  embeddable.  If  B  is  a  single  edge  it  is  easy  to  add  the  edge  to 
the  embedding  of  G'.  Assume  B  is  nontrivial  Since  B  is  nondominating,  v  must 
lie  on  the  exterior  face  /,  of  B  and  ?/,(v)  #  90c  (why?). 

Embed  G'  and  B  separately  and  consider  the  vertex  v  in  both  embeddings 
If  =  180°.  then  v  is  on  only  one  edge  in  G'.  Add-  new  grid  lines  to  the 

embedding  of  G',  create  the  shape  U  as  shown  in  figure  6.2a,  magnify  the  embed¬ 
ding.  and  embed  B  in  the  U  as  in  lemma  5.4a.  If  p/,(v)  =  270°,  then  v  is  either 
on  just  one  edge  in  G‘  ,  or  on  two  perpendicular  edges  in  G'.  In  both  cases,  add 
new  grid  lines,  create  the  shape  V  and  embed  B  as  shown  in  figure  6.2b.  ■ 

Before  we  describe  an  algorithm  for  testing  embeddability,  we  need  an  algo¬ 
rithm  for  testing  whether  "dominate"  is  a  partial  order  on  the  set  of  biconnected 
components  From  the  tree  T  of  biconnected  components,  we  construct  f  a 
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Adding  B  to  the  embedding  of  G‘ 


partially  directed  tree  as  follows  Assume  that  no  biconnected  component  dom¬ 
inates  and  is  dominated  at  the  same  vertex.  If  so  then  "dominate''  is  not  a  par¬ 
tial  order.  Direct  edge  (t i,B)  from  B  to  v  if  B  dominates  at  v.  Direct  edge 
(k  ,B)  from  v  to  B  if  B  is  dominated  at  v .  Leave  all  the  other  edges  undirected. 

This  partially  directed  tree  f  can  be  constructed  in  linear  time  as  follows. 
Find  the  faces  of  each  of  the  biconnected  components  using  the  algorithm 
check-biconnected  This  takes  0(|K|)  time  Check  for  dominations  at  each 
articulation  vertex  as  described  in  definition  6.2.  There  are  at  most  4  bicon¬ 
nected  components  at  each  articulation  vertex  and  hence  there  are  at  most  12 
(ordered)  pairs  to  be  tested  for  domination  (  in  fact  only  2  tests  are  necessary, 
how?  ).  Construct  T  by  directing  the  edges  of  T  as  described  earlier.  Note  that 
articulation  vertices  and  biconnected  components  can  be  found  in  0(\  V|)  [1]. 
For  each  vertex  *  in  T ,  denote  by  <**(*),  dl^(x),  and  d(x).  the  number  of 
incoming  arcs,  the  number  of  outgoing  arcs,  and  the  number  of  undirected 
edges  respectively.  The  rest  of  the  algorithm  is  given  below 

Algorithm  check-domvnate-po(G). 
begin 

construct  T ; 

for  each  vertex  x  in  T  do 

if  «f*,(x)  >  1  then 
begin 

■write  (  Viot  a  partial  order' ), 
quit 

end 

If  search  (f)  then  write  (  <ye s,  partial  order  ) 
rise  write  (  hot  a  partial  order  ' ) 

end 

function  search  (f):  boolean; 

begin 

if  T  =  ♦  then  search  =  true 
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elee  begin 

If  35  £  B  with  (B)=Q,  d+{B)+d(B)= 1  then 

begin 

Let  v  be  the  neighbor  of  B ; 

if  d^t{v)+dovt(v)+d(v)=  1  then  search :  =  search 
else  search  ^  search(T~\B ,t/{) 
end  elae  search  =  false 

end 

The  above  algorithm  can  be  easily  shown  to  be  correct  using  definition  6.3 
and  corollary  6.1.  The  boolean  function  search  can  be  implemented  nonrecur- 
sively  to  run  in  linear  time  by  maintaining  a  queue  of  the  leaves  of  T . 

Given  the  biconnected  components  and  articulation  vertices,  checking  that 
the  articulation  vertices  are  interlace-free  can  be  done  in  0(|K|)  time.  Let 
check-interlace-free  be  a  procedure  that  checks  a  given  articulation  vertex  for 
interlace-freedom.  We  end  this  section  with  a  0(  |  K 1 }  algorithm  for  testing 
embeddability  of  rectilinear  graphs 

Algorithm  check-rectUinear(G): 
begin 

Decompose  G  into  its  biconnected  components, 
for  each  biconnected  component  B  do  check-bicormected(B), 
for  each  articulation  vertex  v  do  check-interlace-free  (v), 
che  ck-dominate-po  ( G) 

end 


7.  An  Embedding  Algorithm 

In  the  previous  section  we  gave  an  algorithm  for  testing  embeddability 
This  algorithm  can  be  easily  modified  into  an  algorithm  which  gives  an  embed¬ 
ding.  However,  the  complexity  of  this  naive  algorithm  would  be  0(  |  F|8).  The  rea¬ 
soning  is  as  follows.  The  path  P'  that  we  find  in  the  proof  of  theorem  5  1  could  be 
0(  j  V ] )  long  For  each  topological  action  that  we  apply  on  this  path  to  transform 
it  to  the  path  P,  we  update  the  coordinates  of  the  vertices  in  the  embedding 
once.  Thus  for  each  path  added  we  require  0{  \  V\2)  time  There  can  be  0(  I  Fi ) 
such  paths  and  hence  the  complexity  of  the  algorithm  is  0(  ]  F|s).  To  reduce  the 
complexity  to  0{  \  V|)2),  we  have  to  make  sure  that  the  path  P'  is  never  longer 
(asymptotically)  than  the  path  P.  In  this  case  the  sum  of  the  lengths  of  all  such 
paths  P'  is  0(  |  K|),  and  the  0(  |  Kj8)  complexity  follows.  In  the  following,  we  show 
how  we  can  always  find  such  paths,  describe  the  algorithm,  and  analyze  its  com¬ 
plexity. 

Lemma  7.1:  Let  0  be  a  planar  biconnected  multigraph  with  minimum  degree 
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three.  Then  any  embedding  of  G  has  an  Interior  face  of  size  at  most  five 
(  Proof:  The  dual  G*  of  G  is  also  a  planar  graph.  Since  G  has  minimum  degree  3. 

G*  is  a  simple  graph.  Hence  G*  has  at  least  two  vertices  of  degree  «  5  [2].  G  is 
biconnected  and  hence  one  of  the  vertices  must  correspond  to  a  face  of  size  £  5 

■ 

I  Lemma  7.2:  Given  an  embedding  of  a  planar  biconnected  graph  G,  which  is  not  a 

cycle,  there  is  a  simple  path  P,  such  that  (i)  the  interior  vertices  of  P  all  have 
degree  2.  (ti)  the  end  vertices  of  P  have  degree  as  2,  (in)  P  appears  in  an  inte¬ 
rior  face  /  in  the  planar  embedding,  and  (iv)  5.  \P  i  \f  | . 

Proof:  As  in  the  proof  of  lemma  5.2,  transform  G  to  G'  by  replacing  all  paths 
-  with  property  (i)  and  (ii)  by  edges.  By  lemma  7.1,  G'  has  an  interior  face  /  of 

size  at  most  5.  The  longest  of  all  the  paths  in  G  corresponding  to  the  edges  of  / 
will  satisfy  conditions  (in)  and  (tu).  ■ 

To  get  an  embedding  of  a  given  rectilinear  graph  we  first  test  if  the  graph  is 
embeddable  and  then  apply  the  following  algorithm. 

I 

Algorithm  embed-rectUinear(G), 
begin 

for  each  biconnected  component  B  do  embed-biconnected  (B). 
join-the-embeddings, 

I  end 

Algorithm  embed-biconnected (B), 
begin 

get-iong-path  (  P,  P,,  a  ); 

I  embed-biconnected  (  B  -  P  ), 

find-path-in-embeddxng  (  P\  Px,  o  ); 
apply-actions-and-transform  (  P,  P  ) 

end 

Procedure  get-long-path  returns  paths  P  ,P\,  and  square  o,  such  that  P 
*  satisfies  the  conditions  of  lemma  7.2,  and  the  interior  face  /  =  PP,  simplifies  to 

a  By  lemma  7.2  such  a  path  exists. 

Procedure  find-path-in-embedding  traces  a  path  P'  in  the  embedding  of 
B-P,  such  that  P'  starts  and  ends  in  the  same  directions  as  P,  and  PP, 
simplifies  to  <r.  P'  and  P  are  both  complements  of  P,  with  respect  to  the  square 
a.  Since  PP,  is  an  interior  face,  P'  can  be  obtained  by  starting  in  the  required 
direction,  then  following  the  path  Pj  in  the  embedding  of  B-P,  and  ending  in 
the  required  direction  (figure  7.1).  This  will  result  in  P'  being  a  complement  of 
Pi  with  respect  to  a.  We  have  \P'  |  =  0(jP,|)  =  G(|P|). 

I  Procedure  apply-actions-and-transform  applies  a  sequence  of  the  four 

topological  actions  to  P'  in  the  embedding  of  B-P+P  and  transforms  it  to  P 
thus  resulting  in  a  embedding  of  B .  This  is  done  by  first  simplifying  the  path  P' 


Figure  7.1 

Finding  the  path  P'  in  the  embedding  of  B—P 
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Figure  7.2 

Path  addition,  simplification  and  expansion 

and  then  expanding  the  simplified  path  to  get  the  path  P  (figure  7.2).  The 
number  of  actions  applied  will  be  0(  | P |  + 1 P’ | )  =  0(\P\) 

Procedure  join-embeddings  takes  the  embeddings  of  the  biconnected  com¬ 
ponents  and  puts  them  together  to  get  an  embedding  for  C.  This  is  done  essen¬ 
tially  following  the  proof  of  theorem  6.1.  Find  a  nondominating  component  B. 

I  Recursively  embed  G'  =  G  -  B.  Join  the  embeddings  of  B  and  G'  using  the 

shapes  U  or  W  as  shown  in  figure  8.2 

The  algorithm  can  be  shown  to  be  correct  using  the  material  developed  in 
the  previous  three  sections,  We  now  analyze  the  complexity  of  each  step  in  the 
algorithm  and  show  that  the  total  complexity  is  0{  j  F|z). 

I  Procedure  join-embeddings  updates  each  coordinate  at  most  once  per 

recursive  call.  The  total  number  of  calls  is  bounded  by  the  number  of  bicon¬ 
nected  components.  Hence  this  procedure  takes  0(  |  K|*)  time. 

Procedure  get-long-path  can  be  implemented  to  run  in  0(|  K|)  time  each 
time  it  is  called  Remember  that  we  can  get  the  faces  of  a  biconnected  graph 

1  from  the  testing  algorithm,  and  searching  all  faces  to  get  the  required  face 

takes  linear  time  Procedure  find-path-vn-emb ■  dding  takes  0{\P\\)~  0(|  V|) 
time.  These  two  procedures  will  be  invoked  at  most  0(|Vj)  time.  Hence  total 


i 


-21  - 


time  spent  in  these  calls  is  0(  |  F|*). 

Procedure  apply-ac tions-and-tronsform.  applies  a  sequence  of  0(\P\) 
actions.  Each  edge  in  G  will  appear  in  only  one  such  path  P.  Hence  the  sum  of 
the  lengths  of  all  such  paths  P  is  0(  |  V)).  Each  action  updates  at  most  0(  |  V|) 
coordinates.  Therefore  the  time  spent  in  calls  to  this  procedure  is  0{  \  F|8). 

B.  Consistent  Rectilinear  Graphs 

Certain  rectilinear  graphs  cannot  be  drawn  on  the  grid  even  if  we  relax  the 
planarity  criterion.  We  say  that  a  rectilinear  graph  G(V,E,X)  is  consistent  if  it 
can  be  drawn  on  the  grid  satisfying  the  ordering  relation  X.  In  other  words.  G  is 
consistent  if  the  set  of  equality  and  inequality  constraints  generated  in  part  1  of 
definition  2.2  is  consistent. 

The  equality  constraints  define  an  equivalence  relation  on  the  set  of  coordi¬ 
nates  of  the  vertices  of  G  Let  us  denote  by  #(x)  the  equivalence  class  contain¬ 
ing  the  coordinate  *.  Denote  by  IM  and  7y  the  sets  of  x-coordinate  and  y- 
coordinate  inequality  constraints  respectively.  Construct  two  directed  graphs 
£%(!*.£*)  and  Gy(Vy.Ey)  as  follows: 

Vi  =  {e(x)  |  x  =  x( a),  a  G  Fj  and  Em  =  \(xlrx2)  J  x,  >  x2  e  /,j 
Vv  and  Ey  are  similarly  defined. 

It  can  be  easily  shown  that  G  is  consistent  if  and  only  if  the  two  directed 
graphs  Gk  and  Gy  are  both  acyclic.  A  solution  to  the  coordinates  which  satisfies 
the  constraints  will  correspond  to  a  nonplanar  embedding  of  G  on  the  grid.  This 
can  be  obtained  by  performing  the  topological  sort  operation  [4]  on  the  two  acy¬ 
clic  digraphs.  In  fact  this  will  yield  a  solution  that  minimizes  the  area  of  the  rec¬ 
tangle  bounding  the  embedding. 

In  a  nonplanar  embedding  of  a  consistent  rectilinear  graph  on  the  grid,  all 
crossings  are  between  horizontal  edges  and  vertical  edges.  The  vertical  edges 
can  be  assigned  one  layer,  and  the  horizontal  edges  can  be  assigned  a  second 
layer.  In  other  words  the  'thickness'  [2]  of  a  consistent  rectilinear  graph  is  less 
than  or  equal  to  two.  A  generalization  of  the  rectilinear  graph  embedding  is  the 
problem  of  embedding  a  rectilinear  graph  with  layers  preassigned,  in  which  no 
two  edges  belonging  to  two  different  layers  cross.  This  remains  an  open  prob¬ 
lem 

0.  Extension*  and  Open  Problems 

As  mentioned  in  the  previous  section,  the  embedding  problem  for  layer 
assigned  rectilinear  graphs  is  still  to  be  solved.  Ibis  has  important  applications 
in  VLSI  It  is  easy  to  show  that  if  a  rectilinear  graph  is  allowed  to  be  discon¬ 
nected,  then  the  optimal  area  embedding  problem  is  NP-complete  (  reduction 
from  two  dimensional  bin  packing  ).  However,  the  question  is  open  for  con¬ 
nected  rectilinear  graphs. 
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On  a  more  theoretical  side,  we  can  define  'triangular',  'hexagonal'  and  other 
polygonal  graphs,  and  consider  the  embedding  problems  on  appropriate  grids. 
However,  it  seems  that  the  work  in  this  paper  does  not  generalize  easily.  This  is 
mainly  because  these  polygonal  graphs  lack  the  nice  simplification  properties  of 
rectilinear  graphs.  Geometry  seems  to  dominate  over  topology  in  these  polygo¬ 
nal  graphs. 

We  conclude  with  a  note  that  our  0(|  VI8)  embedding  algorithm  will  be 
Implemented  in  AL1.  In  ALI,  layouts  are  described  in  a  hierarchical  fashion,  and 
hence  the  algorithm  can  be  applied  hierarchically  on  a  cell  by  cell  basis  In  this 
case  the  complexity  of  obtaining  the  embedding  is  the  sum  of  the  quadratic 
complexity  over  all  cell  instances. 
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Abstract  We  auume  that  long  wire*  represent  large 
capacitive  load*,  and  investigate  the  effect  on  the  area  of 
a  VLSI  layout  when  drivers  are  introduced  along  many 
long  wires  in  the  layout.  We  present  a  layout  for  which 
the  introduction  of  drivers  along  long  wires  squares  the 
area  of  the  layout;  we  show,  however,  that  the  increase  in 
area  is  never  greater  than  this,  if  the  driver  can  be  laid 
out  in  a  square  region.  We  also  show  an  area-time  trade¬ 
off  for  a  single  long  wire  by  which  we  can  reduce  the  area 
of  its  driver  to  9(l,).g<  1.  from  9(1).  if  we  can  tolerate  a 
delay  of  9 (l,_f)  rather  than  9(log  1);  and  we  obtain  tight 
bounds  on  the  worst-case  area  increase  in  general  lay¬ 
outs  having  these  drivers,  using  the  Brouwer  fixed-point 
theorem.  We  also  derive  results  for  the  case  when  drivers 
are  embedded  in  rectangles  that  are  not  square.  Finally, 
we  extend  the  use  of  our  upper-bound  technique  to  other 
layout  problems. 


1.  Introduction 

The  presence  of  long  wires  in  a  VLSI  layout  slows 
down  the  performance  of  the  circuit  due  to  signal  propa¬ 
gation  delays.  This  effect  can  be  minimised  by  using  a 
driver  at  the  head  of  each  long  wire.  However,  the 
drivers  themselves  occupy  space,  and  the  question 
arises  whether  drivers  can  be  introduced  efficiently 
when  there  ere  manv  lone  wires  In  a  lavout.  For  exam¬ 
ple.  consider  the  simple  case  of  a  long  bidirectional  wire. 
A  driver  needs  to  be  introduced  at  each  end  of  the  wire, 
and  the  presence  of  each  driver  increases  the  length 
along  which  the  other  driver  has  to  drive  the  signal.  A 
large  network  with  many  long  wires  may  contain  several 
instances  of  mutual  interactions  of  this  and  other  kinds, 
and  it  is  not  clear  whether  drivers  can  be  Introduced  in 
an  area-efficient  manner  under  these  conditions  (or  even 
whether  they  can  be  introduced  at  all).  This  is  the  prob¬ 
lem  we  analyze  in  this  paper. 

The  question  of  delays  along  long  wires  is  a  complex 
one,  and  no  consensus  has  been  reached  as  yet  on  the 
correct  way  of  modeling  this.  Some  papers  (e.g.,  [BKBO], 
[VBO])  assume  that  a  constant  propagation  delay  models 
the  current  situation  quite  well.  But  it  is  generally 
accepted  that  this  is  not  a  good  approximation  of  the 
physical  situation.  Head  and  Conway  [MCBO]  and 
Chazelle  and  Monier  [CMBl]  suggest  that  a  wire  is  basi¬ 
cally  a  distributed  RC-network,  and  the  delay  thus  is  pro¬ 
portional  to  the  square  of  its  length.  This  delay  can  be 
reduced  to  a  linear  delay  by  using  repeaters  along  the 
wire.  Linear  delay  is  also  the  asymptotic  limit  imposed 
by  the  speed  of  light.  Thompson  (TBO,  761]  suggests  that 
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the  resistance  of  a  long  wire  can  be  ignored  and  that  the 
wire  can  be  viewed  as  a  large  capacitive  load  whose  capa¬ 
citance  is  proportional  to  its  length.  This  view  is  also 
presented  in  [HR79]  and  is  supported  by  Bilardi.  Pracchi. 
and  Preparata  [BPPBl]  who  solved  the  diffusion  equation 
for  a  long  wire  and  suggest  that,  for  predicted  future 
gains  in  technology,  the  wire  can.  indeed,  be  modeled  as 
a  purely  capacitive  load.  For  such  a  load,  it  is  well-known 
that  the  delay  can  be  reduced  to  the  logarithm  of  the 
capacitance  by  introducing  a  sequence  of  drivers  which 
occupies  area  proportional  to  the  capacitance  [MC8G]. 

In  this  paper,  we  assume  that  the  wire  can  be 
modeled  as  a  purely  capacitive  load.  We  address  the  fol¬ 
lowing  question  that  arises  naturally  in  this  context,  but 
has  not  been  examined  before:  What  is  tht  on  the 

area  p f  a  layout  whan  drivers  are  introduced  to  speed  up 
signals  along  many  tong  wires?  We  justify  our 
capacitive-model  assumption  by  appealing  to  the  simula¬ 
tion  results  of  [BPPBl  ].  Further,  if  a  long  wire  must, 
indeed,  be  modeled  as  an  RC-network.  then  the  introduc¬ 
tion  of  repeaters  to  reduce  the  delay  will  increase  the 
area  only  by  a  constant  factor.  This  is  because  such 
repeaters  are  of  constant  size,  and  a  long  wire  may  thus 
be  modeled  as  s  sequence  of  short  wires  connected  by 
nodes  of  fixed  size.  In  the  capacitive  model,  on  the  other 
band,  the  area  occupied  by  a  driver  increases  with  the 
length  of  the  wire  it  drives,  and  thus,  drivers  cannot  be 
abstracted  as  constant-size  nodes.  In  fact,  we  show  that 
there  can  be  definite  area  penalties  when  there  are 
many  drivers  in  a  layout.  Our  results  are  of  particular 
significance  for  upper-bounds  in  area-time  products  for 
VLSI  layouts,  since  most  previously  derived  bounds  have 
ignored  the  delay  along  long  wires.  Either  linear  delay 
should  be  assumed  along  connecting  wires  (this  would 
represent  either  "RC"  wires  with  repeaters,  or  "capaci¬ 
tive"  wires  without  drivers),  or  the  area  expansion 
caused  by  drivers  for  capacitive  wires  should  be  taken 
into  consideration;  alternatively,  some  intermediate 
design  for  drivers  can  be  used  from  the  spectrum  of 
designs  we  suggest  in  section  4.  But  upper-bounds 
derived  using  constant  delay  along  all  wires,  and  no  area 
expansion  for  drivers  do  not  model  the  physical  situation 
well. 

For  a  wire  of  length  I,  the  most  familiar  type  of 
driver  occupies  9(1)  area  while  reducing  the  delay  from 
9(1)  to  9(logt).  In  this  paper,  we  present  a  layout  for 
which  the  introduction  of  such  drivers  along  long  wires 
•guars*  the  area  of  the  layout;  we  show,  however,  that 
the  increase  in  area  is  never  greater  than  this,  if  the 
driver  can  be  laid  out  in  a  square  region.  We  prove  the 
upper-bound  by  a  new  technique  that  uses  a  fixed-point 
theorem;  we  believe  tbs  proof  technique  is  important  in 
its  own  right.  We  also  show  an  area-time  trade-off  for  a 
single  long  wire  by  which  we  can  reduce  the  area  of  its 
driver  to  9(1*  ),g  <1.  from  9(1),  if  we  can  tolerate  a  delay 


of  8(1'"*)  rather  than  8{log  ();  and  we  obtain  bound*  on 
the  woril-cate  area  increase  in  general  layout*  having 
these  drivers,  again  using  the  Brouwer  fixed-point 
theorem.  We  also  examine  the  ease  when  drivers  cannot 
be  laid  out  efficiently  in  square  regions. 

The  paper  is  organizsd  as  follows:  Section  2  presents 
the  VLSI  model  we  use.  In  section  3  we  present  a  graph 
on  n*  vertices  with  6(n )  long  wires  for  which  the  layout 
area  expands  from  6(«*)  to  8(n4)  after  introducing 
drivers  fcr  the  long  wires.  Section  4  examines  the  area¬ 
time  trade-off  for  drivers,  and  sections  5  and  6  use  the 
Brouwer  fixed-point  theorem  to  find  general  tight  upper- 
bounds  for  the  area  expansion  caused  by  drivers  in  arbi¬ 
trary  layouts.  Section  7  generalises  the  technique  to 
other  layout  problems,  and  section  B  concludes  the 
paper  with  a  review  and  some  open  questions. 


2.  The  Model 

We  assume  a  layout  model  similar  to  the  one  used 
by  Thompson[TBO]  and  Leiserson[L90],  The  VLSI  circuit  is 
abstracted  as  a  graph  having  bounded  vertex  degree  and 
it  is  embedded  on  a  planar  grid  subject  to  the  following 
constraints: 

1)  Eech  vertex  represents  a  processing  element  and 
occupies  a  constant  area.  Distinct  vertices  of  the  graph 
are  embedded  at  distinct  grid  intersection  points. 

Z)  Edges  have  unit  width  and  are  routed  along  grid  lines 
with  the  restriction  that  no  two  edges  touch  one  another 
except  possibly  when  crossing  perpendicular  to  each 
other.  Also,  an  edge  cannot  be  routed  over  a  vertex  it 
does  not  connect.  Each  edge  represents  a  connecting 
wire  We  will  refer  to  an  edge  as  a  wire  or  a  line  . 

We  assume  that  we  are  given  such  a  layout  for  a  cir¬ 
cuit  with  certain  edges  identified  as  long  wires.  Note  that 
the  number  and  positions  of  long  wires  is  layout- 
dependent,  and  the  aame  circuit  may  have  another  lay¬ 
out  with  shorter  lines.  We  shall  not  go  into  the  question 
of  derigning  layouts  to  minimize  the  presence  of  long 
lines.  Some  work  has  been  done  on  minimizing  the  length 
of  the  longest  line  in  a  layout  for  certain  classes  of 
graphs  ([PRSBl],  [RSBl],  [BLB2]),  but.  as  we  mention  in 
section  6.  the  minimization  criterion  in  our  case  is  a 
different  one.  which  has  not  been  studied  so  far.  We  will, 
therefore,  assume  that  the  layout  is  given,  and  we  will 
introduce  drivers  where  needed  by  making  local  expan¬ 
sions  without  distorting  the  layout  configuration. 

A  long  wire  is  a  bidirectional  element  electrically, 
but  in  a  VLSI  circuit,  it  usually  connects  active  unidirec¬ 
tional  devices,  so  that  the  signal  always  originates  at  one 
fixed  end  of  the  wire  (which  we  will  call  the  Asad  of  the 
wire)  and  propagates  towards  the  other  end.  Tor  the  rest 
of  thi:  paper,  we  assume  that  wires  come  srith  the  direc¬ 
tion  of  signal  flow  specified,  and  that  bidirectional  wires 
are  deromposed  (conceptually)  as  two  adges  connecting 
the  same  vortices,  but  having  opposite  signal  Bow. 

As  mentioned  in  the  introduction,  we  assume  that 
the  line  Is  purely  capacitive,  and  the  delay  thus  grows 
linearly  with  1,  the  length  of  the  line.  This  delay  ean  be 
reduced  by  using  drivers  to  speed  up  the  signal.  Many 
different  types  of  driver  designs  are  possible,  and  they 
differ  in  the  signal  speed-up  they  offer,  and  the  area  they 
occupy.  We  analyze  these  time-space  trade-offs  in  sec¬ 
tion  4.  For  the  present,  we  assume  the  most  familiar 
type  of  driver,  which  reduces  the  delay  to  B(log  I )  at  the 
expense  of  occupying  •(!)  area.  We  also  assume  that  the 
driver  can  be  laid  out  in  a  square  of  olds  k  VT,  for  some 


figure  1 

A  driver  introduced  along 
part  of  a  long  wire 

constant  k.  This  is  not  strictly  true,  since  the  driver  has 
a  transistor  whose  channel  is  B(t )  wide,  and  hence  it  is 
more  naturally  laid  out  in  a  long,  narrow  rectangle. 
Clever  design  methods,  however,  such  as  the  use  of  a 
zig-zag  poly  line  tor  the  gate  in  MOS  technology,  can  be 
used  to  overcome  this  problem.  We  have  more  to  say  on 
this  point  in  section  6. 


figure  2 

Introducing  a  driver  by  making 
cuts  in  a  layout 

Given  a  long  line  w  of  length  I  in  a  layout,  an  area- 
efficient  way  to  achieve  signal  speed-up  would  be  to 
introduce  a  long  skinny  driver  (the  natural  design)  along 
a  portion  of  w  (see  figure  1).  However,  VLSI  layouts  are 
dense  in  general,  and  the  above  scheme  will  not  work  if 
there  are  vertical  line*  present  in  this  region.  In  such  a 
case,  we  make  two  "cuts"  or  "slices"  (a  notion  introduced 
by  Leiserson  [100])  at  the  heed  A  of  each  *  vT  wide, 
one  each  in  the  vertical  and  horizontal  direction*  (figure 
2).  Any  edge  that  is  broken  by  the  cut  is  then  joined 
across  the  cut  by  a  straight  tins.  This  forms  a  square 
region  of  k  vT  side  at  A  with  no  edges  or  vertices  from 
the  original  layout.  Wt  introduce  the  driver  in  this 
region.  This  construction  does  not  disturb  the  layout 
configuration,  and  introduces  the  driver  by  local  expan¬ 
sions.  It  increases  the  length  of  each  side  of  the  layout 
by  the  size  of  the  cut  A  also  increases  by  k~Sl .  (he 
long! A  of  eacA  horizontal  wire  in  f  As  layout  which  inter¬ 
sects  IAs  vertical  tins  drawn  through  A  and  sac  A  vsrftraf 
wirs  that  Wilerescts  IAs  Aorise nlaf  I# is  drown  through  h 
These  other  wire*  will  now  require  larger  drivers  than  the 
ones  they  would  have  needed  in  the  absence  of  a  driver 
at  A.  The  introduction  of  drivers  at  these  other  tines  will. 
In  turn,  affect  the  lengths  of  more  lines,  and  could  even 
increase  the  length  of  w  so  that  the  area  of  Its  driver 
would  have  to  be  revised  upward. 
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We  are  now  in  •  position  to  state  the  problem  we  are 
going  to  analyze: 

ftwn  a  VLSI  layout  with  cerium  edges  specified  as 
long  wires.  what  is  tha  increase  in  the  area  of  tha  layout 
in  tht  wont  ease  whan  drivers  are  used  to  speed  up  sip* 
note  along  thata  long  wires? 

3.  A  Worst -caae  Example 

Consider  the  graph  layout  of  figure  3.  It  consists  of 
a  square  mesh  on  n*  vertices  (assume  n  even),  together 
with  n/2  long  horizontal  lines  W|.  wt,  ■  ■  ■  ,u^/t  If  we 
assume  the  vertices  to  be  laid  out  in  an  nan  grid  (there 
will  be  a  slight  expansion  needed  to  accommodate  the 
long  wires  along  grid  lines,  but  we  can  ignore  this  for 
"order-of-magnitude"  arguments)  and  let  m «  n/2,  then 
so,  runs  from  vertex  (m.l)  to  vertex  (n.l),  wt  runs  from 
(m-1  2)  to  (n-1.2).  •  •  •  .  and  tit,  runs  from  (l.m)  to 
(m  4 1  to  ).  The  signal  goes  from  left  to  right  along  each  of 
these  of  these  lines. 


figure  3 

A  graph  layout  with  0(n2)  area  and  0{n)  long  lines 
of  length  0(n)  each 


We  consider  now  the  increase  in  the  area  of  the  lay¬ 
out  when  we  introduce  drivers  on  wires  U|  to  wm.  Notice 
that  the  presence  of  the  underlying  square-mesh  struc¬ 
ture  forces  us  to  make  cuts  to  introduce  the  drivers. 
Since  n/2  drivers  are  introduced  diagonally,  and  the 
length  of  each  wire  is  0(n)  to  start  with,  the  length  of 
each  ride  of  the  layout  increases  to  at  least  0(n  n'/s),  so 
that  the  area  of  the  layout  increases  from  0(n*)  to  at 
least  0(n’).  However,  this  does  not  account  for  the 
increase  in  the  lengths  of  some  wires  due  to  the  intro¬ 
duction  of  drivers  on  other  wires.  We  now  derive  a  non¬ 
linear  recurrence  relation  that  accounts  for  this  effect, 
and  use  it  to  bound  the  increase  in  area  as  0(n4). 

Since  the  signal  travels  from  left  to  right  along  each 
of  the  long  wires,  the  length  of  W)  is  not  affected  by  the 
introduction  of  drivers  at  the  beads  of  the  other  long 
wires.  In  general,  the  length  of  w,  increases  by  exactly 
the  lengths  of  the  drivers  introduced  on  «/,  to  If  we 
let  d,  be  the  width  of  each  side  of  the  driver  tor  wire  u>(, 
then  we  have 

d^km1'* 
d,«t  (m+d,),/* 


/•i 


Let  s,  be  a  new  variable  defined  by 
e,*» nv-d, 

Sjvm+dj+d, 


at  *m.+  £d, 
f»i 

>a(.,+4 

m«i-i+*(m*£dj) 


i/z 


Tor  convenience,  define  itsm.  Note  that  s,+i  gives  the 
width  of  the  layout  from  the  rightmost  end  to  the  left 
and  of  the  ith  driver,  and  sn  +m  gives  the  length  of  each 
aide  of  the  layout  after  the  introduction  of  drivers  on  all 
of  the  long  wires.  We  find  bounds  for  this  value  below: 


Lemma  There  exist  constants  kt  and  kt  such  that  for 
any  integer  m  >0. 

kt<*+m  ks,k  k,i*+m  for  t  »'/m. 


Proof  We  prove  both  aides  of  the  inequality  by  induction. 
Assume  s,tk  ,t*+m 

Since  s0sm.  the  assumption  holds  for  sc  lor  all  values  of 
A).  Assume  the  result  holds  up  to  s(.  Then. 


fck  ,i»+m  +k  (fc  ,t*+m  )w * 

*k ,  (i  + 1  )*♦  m  ♦  (fck }'  *  -2k ,  )i  -k , 


Hence.  e,«)»m+fc,(t4.1)1  for  k,<(k*/4)  and  for 
sufficiently  large  i.  For  k  ,*(**/  B),  the  result  holds  for 
all  values  of  t.  So  we  have  s,t(t*/  9)i*+m  for  all  ikO. 

Assume  siakti*+m 

Assume  the  result  holds  up  to  a,.  Then, 


*k  (k 

*k tt*+m-rfc((k,vl)s*)1/*  for  ifcm 1/2 

**»(t  +  l)*+m  -(aiI-k(*t»l)l/i<H^i 

<kt(i  +  l)2+m  for  kgfc(k2/B)  (lv^/^♦167fc,). 


Hence  we  have  a1*(fc*/B)(l  +  Via  16/ k*)i*+m  for 
tav'm  .  For  i«Vm  ,  we  need  ty/^*ktm*m,  and  since 
SJZ:  <  niiiVw  Vm  ♦  2kjl  *k)m  .  the  inequality  holds 
when  *(fcfc  Vj  +  2*+z*«. 

Therefore.  s,*k,i,*m  for  all  ikv'rn  when 
*gkmaz(  (k */ 6)( V l ♦  1 6/T*) .  k  Vlv2k  *2kw  ) 

The  required  result  follows  from  combining  the  two 
inequalities.* 

Thus,  the  width  of  the  layout  after  the  introduction 
of  drivers,  which  is  given  by  sm*m,  is  0(m2)=9(n2)  end 
the  area  thus  grows  to  exactly  0(n4).  the  square  of  the 
initial  area.  Tbe  width  of  the  driver  tor  u,  is  given  by: 


ekv'i;'" 

■•(\ ).  tor  t»Vm  . 


1 


Hence  there  are  0(n)  drivers  of  width  0(n)  in  this  lay¬ 
out.  Note  also  that  many  of  the  short  wires  have  their 


lengths  greatly  increased  by  the  introduction  of  drivers 
for  wires  uit  to  tu*.  For  instance,  the  wire  connecting  the 
vertex  (l.n)  to  the  vertex  (2,n )  has  its  length  increased 
from  1  grid  unit  to  9{n)  grid  units.  We  may  now  want  to 
introduce  drivers  for  thesr  wires,  too.  but  since  we  a 
deriving  a  lower  bound  in  this  section,  we  will  ignore  this 
point. 

The  layout  of  figure  3  is  simple  in  the  sense  that  it 
has  no  mutual  interaction:  there  do  not  exist  wires  uy 
and  tu,  such  that  the  driver  of  tu,  increases  the  length  of 
tu,  and  the  driver  of  increases  the  length  of  su,.  In 
general,  we  could  have  a  more  complex  pattern  of  long 
wires  in  which  such  mutual  interactions  occur  in  many 
ways.  However,  in  section  6.  we  present  the  surprising 
result  that  the  penalty  in  area  never  becomes  worse 
than  the  square  of  the  original  area. 

4.  Area-Delay  Trade-ofl  for  a  Driver  of  a  Long  Wire 

When  a  large  capacitive  load  cL  is  driven  by  an 
inverter  with  resistance  r,.  the  delay  is  the  product  r4c i 
(note  that  r,  here  refers  to  an  "average"  resistance, 
since  the  resistance  and  the  delay  vary  depending  on 
whether  the  input  is  low  or  high.  In  what  follows,  we 
make  some  other  simplifications  of  the  physical  situa¬ 
tion,  but  these  will  affect  the  constant  factor  only,  ao 
that  tne  "order-of-magnitude"  results  bold).  If  the  resis¬ 
tance  and  capacitance  of  a  minimum-size  inverter  are  r, 
and  c,  respectively,  then  the  basic  unit  of  delay  is  l«r,c, 
and  this  represents  the  delay  when  a  minimum-size 
inverter  drives  the  gate  of  another  minimum-size 
inverter.  In  terms  of  this  unit  t ,  a  minimum-size  inverter 
driving  a  load  capacitance  ct  causes  a  delay  of  cLt/c,. 
and  this  grows  linearly  with  the  load  capacitance.  When 
the  load  is  a  long  wire,  this  delay  grows  linearly  with  1, 
the  length  of  the  wire,  assuming  the  width  and  thickness 
of  the  wire  remain  fixed. 


- g>-~ 

e  i*/  i'e,  es=/  r/sei 

ri=r,//,  *-t=r,//,/t 


channel  width  of  inv.  1  *  /  ,•  (minimum  channel  width) 
channel  width  of  inv.  2  s  /,  (channel  width  of  inv.  1) 

figure  4 

Delay  on  one  stage  of  a  driver  sequence 


To  reduce  the  delay  to  more  acceptable  levels,  we 
use  e  sequence  of  drivers  with  successively  wider  chan¬ 
nels  to  drive  the  load  capacitance.  To  understand  why 
this  scheme  works,  consider  the  delay  when  one  inverter 
with  channel  width  f  \  times  the  minimum-size  inverter 
drives  the  gate  of  another  inverter  whose  channel  is  /s 
times  wider  than  that  of  the  driving  inverter  (see  figure 
4).  The  resistance  of  the  first  inverter  is  r,//t  and  the 
capacitance  of  the  second  inverter  is  /i/Ki  end  hence 
the  delay  of  this  stage  is  /  ft  .  If  we  now  have  n  Inverter 
stages,  starting  with  a  minimum-size  inverter  as  the  first 
stage,  a  driver  with  a  very  wide  channel  as  the  nth  stage 
to  drive  et,  and  intermediate  stages  with  successively 
wider  channels,  then  the  total  delay  is  the  sum  of  the 
delay  through  the  n  stages.  Let  f,  be  the  ratio  of  the 
channel  width  of  the  t+lst  inverter  to  that  of  the  itb 


inverter  for  i“  1,2,  -  -  ■ .«— 1,  and  let  /„  be  the  ratio  of  c i 
to  the  capacitance  c.  of  the  last  driver  (figure  5). 


channel  width 
(normalized)  1  / , 


figure  5 

A  sequence  of  drivers  for  load  cL 

Ve  have  /,/,  -  -  -  /,*ej/C|  and  the  total  delay  through 
the  n  stages  is  T*(f  \*fi*  •  ■  ■  ♦/«)!•  If  we  fix  each  /, 
to  be  equal  to  a  constant  /,  then  the  delay 
T*log/(c|,/c1),  and  this  is  the  minimum  delay  we  can 
obtain  up  to  a  constant  (see  [kCBO]).  However,  the  area 
occupied  by  the  drivers  in  this  case  grows  linearly  with 
eg,  as  we  show  below,  and  it  may  be  desirable  to  reduce 
this  area  penalty  at  the  expense  of  a  slightly  greater 
delay.  Ve  show  below  that  such  a  space-delay  trade-off 
exists.  Some  of  these  results  are  implicit  in  the  work  of 
Lin  and  Unholm  [LL75]  and  Jaeger[J75]. 

Let  the  area  occupied  by  a  minimum-size  inverter 
be  A,.  Then,  an  inverter  whose  channel  is  /  times  as 
wide  as  a  minimum-size  inverter  requires  an  area  of 
about  /A,.  Hence  the  total  area  A,  occupied  by  the 
sequence  of  drivers  with  logarithmic  delay  is  given  by 

A*\  (!♦/♦/•♦  4/**,> 

«e  (/*) 

«e(cj.). 

Now,  consider  the  general  sequence  of  drivers  described 
above.  The  area  occupied  in  this  case  is 

*«*(!♦/ ,♦/,/*♦  ♦/,/■  /.-l) 

Assume  that  we  can  tolerate  a  delay  of  9(i*).0<q<l 
(*9((ci)*)).  Since  the  area  of  the  last  driver  dominates 
the  above  expression  for  area,  we  look  for  the  amallest 
size  for  the  last  driver  which  will  still  result  in  a  delay  of 
not  more  than  9(C).  Since  c**^//*.  and  the  delay 
through  the  last  stage  is  /«< .  the  best  we  can  do  is  to  set 
/**9(C),  and  this  results  in  an  ares  of  9(1' "*)  for  tn. 
Now,  the  question  is  whether  we  can  accommodate  the 
other  drivers  to  occupy  not  more  than  9(1'"* )  area  end 
cause  no  greater  than  9(C)  delay.  The  answer  to  this  is 
'yes',  and  in  fact,  we  can  accomplish  this  in  many  ways 
Vs  describe  one  method  of  doing  this,  which  minimizes 
the  constant  on  the  delay  term.  For  this,  we  make  t„ 
the  load  so  that  we  need  to  drive  a  load  of  9(t'**).  Ve 
know  from  the  design  for  drivers  with  logarithmic  delay, 
that  this  can  be  done  w.tb  9((l-f  )logi)  delay  and  area 
9(1,_»)  using  9((iog  r„/c())  stages.  Hence  the  total  arcs 
for  the  n  stages  remains  9(1*'*)  while  the  delay  is  incre¬ 
mented  by  only  a  9(log  1)  term  and  bence  remains  9(C). 

Ve  have  thus  shown  the  following  interesting  area- 
delay  trade-off  for  the  driver  of  a  single  long  wire  of 
length  <  :  For  any  7  0<q  <1,  sue  eon  design  a  driver  urjh 
delay  9(C)  and  are a-detay  product  9(1 )  If,  Aosuevrr.  u.c 
need  to  reduce  (As  delay  to  9(iog!)  suAicA  is  the 
minimum  delay  acAievahit,  (Asn  (As  area-delay  product 
gees  up  to  9(1  log  I ) 
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6.  An  Upper-Bound  Proof  Technique  Using  Fixed-Points 

Let  f(*,.zt.  -  -  -  ,xm)  be  •  mapping  from  IP  to  R“ 
defined  by 

F(x,,  ■  •  '  ‘  •**)*  ‘  ‘  •/■»(*!»  ‘  ‘  •*■•))» 

where  each  /,  ie  a  mapping  from  R*  to  R  Than,  a 
/ixttf  -point  of  f  it  a  rector  (*,.•■  je*  )eR“  tucb  that 

/■(*,.  •  ■  .**)*(*.■  '  '  -O. 

/«(*..  .a.)**,.  for  <*1.  ■  ,m. 

The  following  are  tome  well-known  remits  from  real 
analysis  (tee  e.g..  Bartle  [B84]). 

Brouwer  RiedPoint  Theorem  Let  B>0.  and  let  D  be  the 
hypersphere  defined  by  D  *  t*£R“:|* \*B{.  Then,  any 
continuous  function  with  domain  D  and  range  contained 
in  D  has  a  fixed-point. 

Pact  1  Any  real  continuous  function  on  a  bounded  closed 
eet  in  R*  bas  a  maximum  and  a  minimum. 

Pact  2  For  B>  0,  the  hypersphere  defined  by  D  * 
(z  cR*‘  |*|  t.B\  it  a  closed  bounded  tet. 

be  win  nerd  the  following  two  definitions: 

Definition  )  Let  /(z,.  ,xm)  be  a  continuous  mapping 
from  HP  to  R  and  let  g  (*)  be  the  mapping  from  R*  to  R* 
such  thit  for  each  R*0.  g(ff)  gives  the  maximum  abso¬ 
lute  value  of  /  in  the  closed  hypersphere  of  radius 
l?vn  .  He  call  g  the  behavior  function  of  /  .  Note  that  by- 
facts  1  and  2.  g  always  exists. 

Definition  2  A  function  g  (z )  from  R*  to  R*  il  wvblintar 
in  z  if  there  exists  sn  r*0  such  that,  for  all  **r.  we  have 
g  (z  )iz  We  call  r  a  breakpoint  of  the  function  g . 

We  now  prove 

Theorem  1:  Let  F  be  a  continuous  mapping  from  R“  to 
given  by 

/•(*,.  •*»)*(/  i(*i.  •*».)•  ./•(*».  •*«.»• 

and  for  each  /,,  let  the  behavior  function  g,  be  eublinear 
in  *.  with  breakpoint  r,.  Then, 

1)  the  function  F  has  a  fixed-point  within  a  hypersphere 
of  radius  /fVm  .  where  /?=max  rt,  and 

2)  each  component  of  the  fixed-point  is  no  greater  than 
R  in  magnitude 

Proof  To  prove  1).  it  it  sufficient,  from  Brouwer’s 
theorem,  to  prove  that  the  range  of  F  it  contained  on  or 
within  the  prescribed  hypersphere  when  the  domain  it 
similarly  constrained.  For  this,  it  is  sufficient  to  prove 
that,  for  each  i,  gx(R)*R.  since  this  will  guarantee  that 
the  Euclidean  length  of  the  range  is  no  greater  than 
>/mR  But  this  follows  immediately  from  definition  2, 
since  the  functions  g ,  are  sublinear  with  breakpoints  r«, 
and  for  aach  *.  tt*R. 

Note  that,  by  requiring  g,(R)*R  for  each  i,  we  are 
constraining  the  range  to  tie  within  a  closed  hypercube 
inscribed  within  the  hypersphere  of  tbe  domain.  This 
additional  restriction  ensures  that  each  component  of 
the  fixed-point  ie  at  most  R  in  magnituds,  and  this 
proves  2),  and  tbe  theorem  • 

In  our  application  to  layout  problems,  and  In  partic¬ 
ular,  to  tbe  problem  of  drivers,  we  will  be  interested  in  a 
special  class  of  functions  /,(*,.  -  -  -  .a* )HP-*R  dtflnad 


by 

/«(*„.. wben(f,+ 

/•i  f*i 

■  0,  otherwise.  (1) 

where  k  and  1,  are  positive  constants,  0<p  <  1 .  and 
0 BjyWAf,  for  some  positive  constant  K.  Usually.  4V  it 
either  0  or  1.  For  aucb  functions,  we  have  the  following 
corollary  to  Theorem  1. 

Corollary  1.1:  Let  F  be  a  mapping  from  R“  to  R“  given 
by 

F(*  !,■••*,)•(/  ,(*J.  •*„).  ./*(*!■  .*w)). 

where  the  s  are  given  by  (1)  Then  F  bas  a  fixed-point, 
each  component  of  which  is  positive  and  bas  magnitude 
no  greater  than  R.  where 

*-ma xfk.ff,)*.  .m"'1-**).  the  k.’s 

being  positive  constants  independent  of  m . 

Proof  h  is  easy  to  prove  that  aach  /,  has  behavior  func- 
tion  gt (*)=/,(*,  ur).  Let  r,*(a  solution  of  tbe  equa¬ 
tion  k (It*  £4,,*  )***)•  Then,  since  p<l.  g,  is  sublinear 

/« i 

with  breakpoint  r,.  Hence  from  Theorem  1.  each  com¬ 
ponent  of  a  fixed-point  of  F  baa  magnitude  no  more  than 
B.  where  J?*mex  rt. 

Cast  1  Assume  that  each  l,*(£dMr,)  Then  we  require 
ft 

e,(^  4tt)>(r,V'er,.  tor  aome  conxtant  e,«2 *.  or, 
1“  _ 

r«*c'«(2^v>>/ll",>’  wb*re  Tbe  maximum 

postibh-  value  for  erch  4^,  is  K  Hence,  each  coordinate 
of  the  fixed-point  is  guaranteed  to  have  a  value  no 
greater  then  R  *km  for  aome  constant  km4, 

(that  depend!  on  k  end  K). 

Chat  2  There  exist  f,’s  such  that  aach  1*>(  l  4 *)r» ) 

r*i 

Let  4  be  the  maximum  value  of  such  l*’s  Then 

£  4Vr*“  0(0.  and  benca  we  need  for  some  con- 

j»i 

stant  *,<2 k.  Hence,  in  this  case,  each  coordinate  of  the 
fixed-point  is  guaranteed  to  have  a  value  at  most 
tfama x(k.  If  .*m. 

By  combining  the  two  cases,  the  required  result  fol¬ 
lows.  Clearly,  since  the  range  of  F  lies  entirely  in  the 
positive  orlbanl,  each  component  of  the  fixed-point  must 
be  positive.* 

Corollary  1.2  If  at  most  r  of  the  4W  can  be  nonzero  for 
oacb  function  /,  In  (1),  then  each  component  of  a  fixed- 
point  of  F  is  positive  and  bat  magnitude  no  greater  than 
R,  where  Ramax(fc , if.  ,km  lfc.k«Mr*''<,"')) 

Proof  Immediate  from  the  proof  of  Corollary  1.1.* 

I.  General  Upper  Bounds  for  Area  Penalty  ta  layouts 
with  Drivers 

We  initially  assume  that  a  driver  for  a  wire  of  length 
1  can  be  embedded  In  a  equare  of  aide  kl*  (p*  1/2) 
When  pa>/2,  this  corresponds  to  the  standard  driver 
with  logarithmic  dalay.  Since  we  are  new  proving  upper 
bounds,  we  wit)  assume  all  wires  in  the  layout  to  be  long 
TTiis  may  be  a  necessary  assumption,  since  a  wire  that 


»» 


was  short  Initially  may  become  long  due  to  the  introduc¬ 
tion  or  driven  on  other  wirei  (a>  we  aaw  in  aection  3). 

Now.  aaaume  that  the  layout  la  in  an  nan'  (rid 
(nan  ).  Hence  the  area  of  the  layout  before  introducin( 
driver*  is  A*nn'.  We  note  the  following  points  about  long 
line*  and  Uiair  drivers  in  this  layout: 

1)  If  each  long  line  is  f  units  long,  then  there  are  at  most 
2 A/ 1  long  lines  in  the  layout. 


I 

4* 

-<• - 4* - >- 

i 

t 


figure  6 

The  position  of  a  cut  for  a  driver 

2)  Thr  cut  perpendicular  to  the  direction  of  the  wire  is 
always  made  with  thr  position  of  the  head  fixed,  and  on 
the  side  the  wire  extends  (figure  €j.  Thus,  if  there  are 
two  wires  with  their  heads  along  the  same  vertical  (or 
horizontal)  line  end  with  opposite  signal  flow,  their 
drivers  do  not  oveilap  the  same  horizontal  (or  vertical) 
region .  Since  this  simplification  at  most  doubles  the 
effect  of  drivers  at  any  single  vertical  or  horizontal  posi¬ 
tion.  the  "order-of-magnitude"  results  remain 
unaflected.  The  cut  in  the  direction  parallel  to  the  wire 
is  made  on  the  side  of  the  wire  which  lies  on  a  clockwise 
rotation  from  the  wire  (see  again  figure  6).  With  this 
convention,  we  can  introduce  drivers  on  all  four  lines  at 
a  vertex,  if  needed. 


driver*.  A  similar  result  bolds  for  horizontal  cuts.  Note 
that  this  means  that  we  need  make  at  most  two  cuts  at 
each  vertical  and  horizontal  grid  position  (the  factor  of 
two  arises  from  the  convention  we  adopted  in  observa¬ 
tion  2  above,  regarding  the  position  cf  the  cut).  Thus,  at 
most  2(nvn  )  out*  will  be  needed  to  accommodate 
drivers  in  any  n  xn  layout.  The  widths  of  the  cuts  need 
to  be  determined,  and  we  bound  these  values  later  in 
this  aection  to  obtain  tight  upper-bounds  on  the  worst- 
ease  area  expansion  caused  by  the  introduction  of 
drivers. 

4)  If  a  wire  bends  (i.e..  goes  both  vertical  and  horizontal), 
the  analysis  still  holds.  We  assume,  however,  that  a  wire 
bends  at  most  a  constant  number  of  times,  and  hence  a 
wire  can  have  length  no  more  than  An,  for  aome  con¬ 
stant  K- 

6)  If  a  wire  is  1  grid  units  long,  then  its  length  is  affected 
by  at  most  21  other  drivers.  This  worst  case  happens 
when  there  are  the  heads  of  tiro  long  wires  (either  verti¬ 
cal  or  horizontal)  at  each  of  the  1-1  Inner  grid  positions 
of  hi,  and  the  signal  flow  on  these  wires  is  directed  so 
that  cuts  are  required  on  both  sides  of  the  grid  position: 
and  there  is  a  long  wire  at  each  end  grid  position  of  xl, 
and  aacb  of  these  wires  requires  a  cut  on  the  side 
affecting  the  length  of  wire  in .  This  situation  is  shown  in 
figure  8 


length  of  wire  in  *  3  units 
figure  B 

Worst-case  effect  of  other  drivers  on  the  length  of  o  wire  in 
initially  k  grid  units  long 
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figure  7 

The  effect  of  long  wires  with  their  with  heads  aligned 

3)  Consider  the  situation  in  figure  7.  All  of  the  lines 
itVioj.iut.m,  have  their  beads  along  the  vertical  line 
through  «,| .  and  the  signal  flow  on  each  wire  is  such  that 
a  cut  is  required  immediately  to  the  right  of  in,  to 
accommodate  the  driver  for  each  wire.  In  such  a  case, 
the  width  of  the  widest  cut  Is  sufficient  to  include  all 


Consider  a  general  layout  with  m  long  wires 
in  ,.int.  -  .in*  of  initial  lengths  !|.!|.  .1*.  As  men¬ 

tioned  earlier,  bidirectional  wires  ore  treated  os  two 
separate  wires,  with  opposite  signal  flow,  bet  the  driver 
for  the  <tb  wire  be  embedded  in  a  square  of  side  d,  he 
need  to  solve  tor  the  d,'s,  4*1.2.  ,m.  Tor  this,  we 

note  that  the  length  of  the  ith  wire.  is  increased  to  a 
new  value  f,'  by  the  presence  of  drivers  on  some  of  the 
other  wires,  i.e.,  we  have 

Vw^eJdjjd,,  .m. 

»»i 

where  <v  is  nonzero  if  the  driver  of  wire  in,  increases  the 
length  or  in,,  and  is  0  otherwise.  Since  we  allow  a  wire  to 
bend  a  constant  number  of  times,  a  single  cut  can 
expand  a  wire  by  up  to  A'  times  the  width  of  the  cut.  and 
hence  dv  Is  a  nonnegative  integer  having  value  no  more 
than  K  Now.  the  driver  of  in,  will  increase  the  lengths  of 
all  horizontal  wires  that  intersect  the  vertical  line 
through  the  bead  of  it^.  and  all  vertical  wires  that  mter- 
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■ect  Use  horizontal  line  through  the  head  of  iuj  (tee 
figure  9). 


The  tignal  flow  can  be  in  either  direction  in 
each  of  wires  it..  to  w, . 

figure  9 
Determining  6y 

Hence  a  simple  preprocessing  of  the  layout  gives  the 
values  for  *1.  -  -  ,m.  Note  that  we  can  treat  both 

horizontal  and  vertical  wires  uniformly  because  of  our 
assumption  that  the  driver  occupies  a  square  region. 
Later,  we  will  need  to  modify  our  equations  when  we  con¬ 
sider  drivers  embedded  in  rectangles  that  are  not 
square.  The  value  for  each  d.  it  given  by 

*-*ar 

(2) 

Let  us  rewrite  equation  (2)  by  defining  d,  in  a  func¬ 
tional  form,  and  extending  its  domain  to  R“  as  follows: 

d,  =  Md,.  ■Ox*(I.+  £4vdi)'.  i  =  l,  .to.  when 

<•  i 

=0  otherwise. 

Let  F(d  |.  ■  •  .dm )  he  the  mapping  from  RP  to  defined 
by 

F(dt.  .**)«(/ i(*i.  .dm).  ./mid,.  •  .*«))• 

Then  the  solutions  of  the  system  of  equations  (2) 
correspond  exactly  to  the  fixed  points  of  the  function  F. 
i.e  .  tc  the  values  of  d,.  •  ■  ■  .dm  such  that 

Fid,.  .dm)*(d>.  .dm) 

From  the  results  of  the  previous  section,  we  can  now 
immediately  derive  upper  bounds  on  the  width  of  each 
driver.  We  first  determine  this  for  the  case  when  the  lay¬ 
out  is  in  a  square  nxn  grid.  We  later  show  that  this  is  the 
worst  case. 

Assume  tbet  each  f,  is  (  units  long,  and  hence,  from 
observation  6.  up  to  21  of  the  dv  can  be  nonzero  for  each 
i.  Hence,  by  Corollary  1.2.  each  component  of  the  fixed- 
point  has  magnitude  no  greater  than  F.  where 

*«max(k,l*.  •k.l'.fc.,,^  )>'<"»») 


become  greater  than  cl*/(l*9)  for  some  constant  c.  The 
maximum  width  occurs  in  the  case  when  i*8(n)  (since 
this  is  the  maximum  value  for  1.  by  observation  4).  and 
this  width  is  bounded  by  The  worst-case  over  p 

occurs  when  p*l/2  and  in  this  case  the  width  is  0(n). 
But  we  know  from  the  example  of  section  3  that  a  driver 
can  indeed  become  as  sride  as  8(n).  Hence  this  bound  is 
tight.  There  can  be  up  to  8(n)  wires  in  the  layout  (by 
observation  1).  and  hence  we  may  have  to  make  8(n) 
cuts  in  each  of  the  two  directions,  so  that  the  layout  area 
will  increase  from  A*n*  to  8(n4)*8(A*).  This  upper 
bound  in  the  increase  in  area  is  again  matched  by  the 
lower  bound  of  section  3.  For  p<  1/2.  the  area  increases 
to  0((n‘/<‘-*>)»)=OU,/!,-7>). 

If  l«8(nT),0<r  <1,  then  the  width  of  each  driver  is 
0(n’»/!*-»>),  but  by  observation  1.  there  can  be  up  to 
8(n*"T)  such  lines  in  the  layout  We  thus  have  more  than 
8(n)  drivers  in  this  case,  however,  by  observation  3.  if  we 
znake  two  cuts  of  width  8(nr)  at  each  vertical  and  hor¬ 
izontal  grid  position  where  drivers  are  needed  (there  are 
at  most  4n  such  cuts),  then  we  can  accommodate  all 
drivers.  Each  side  of  the  layout  will  expand  to  at  most 
and  the  area  grows  to 
8(n,*,'7/t,■»,))=8(A,•,7/I,■7,)•  Since  r<l.  this  increase  is 
less  than  the  increase  for  the  case  when  r»l. 

If  we  had  a  mixture  of  wires  of  varying  length,  the 
width  expansion  is  still  by  the  same  argu¬ 

ment.  In  the  worst  case,  we  may  have  to  make  a  cut  for 
a  driver  at  each  side  of  each  vertical  and  horizontal  grid 
position.  There  are  4n  such  vertical  and  horizontal  grid 
positions,  and  the  worst  case  occurs  when  B(n)  cuts  of 
width  are  required  in  each  of  the  two  direc¬ 

tions.  We  have  thus  proved  the  following  theorem: 

Theorem  2  If  the  driver  of  a  long  wire  of  length  1  can  be 
embedded  in  a  square  region  of  side  8((r).0<pcl/ 2. 
then  the  area  A  of  any  VLSI  layout  embedded  on  a  square 
grid  increases  to  at  most  0(A,/ll'7))  with  the  introduc¬ 
tion  of  such  drivers  along  long  wires. 

We  have  already  shown  that  this  bound  is  tight  when 
p- 1/2.  We  now  show  that  this  holds  for  any  pel/ 2. 

Lemma  If  the  driver  of  a  long  wire  of  length  1  can  be 
embedded  in  a  square  region  of  side  8(!*).0<p*l/ 2. 
then  there  exists  a  VLSI  layout  whose  area  increases 
from  A  to  n(Al'l,'7))  when  such  drivers  are  introduced 
along  long  wires 

Proof  We  use,  once  again,  the  layout  of  figure  3.  The 
recurrence  relation  for  *,  is  now 

with  se*mvn/2.  As  before.  am+m  gives  the  width  of 
each  side  of  the  Isyout  after  the  introduction  of  drivers 
We  prove  by  induction  that  sm  *  n(TO,/l'"*)). 

Assume  s,  a  e  t  1/<,'7>4-to. 

The  result  bolds  for  So  for  all  values  of  c.  Assume  that 
the  result  bolds  up  to  s,.,.  Then. 

a,  «*,-,♦*  {*,_,>*’ 

«c(<-l)1/l,-*4m+k<c(i-l)>/<'-»>*m)' 
we  (i -1),/I'-7>vto  ♦ke>{i-l)"<'-7' 

■  ci‘/f'*»»((l-l/i)i''“->> 

♦  (k/(c*-7i))(l-l/<)"',-*>)+m 


Hence  the  sridth  of  any  driver  In  this  layout  will  never 
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Since  pk  1/2,  we  have  p/ (l-p)*l,  and  hence 
Hence, 

s.*ct1'<1'*>(l-l/i)(H-((*/e'-*)-l)/<)-i-m 
fceii/(|-»)+m  vbenihZ  end  c*{k  / 

We  elio  need  *j=*n+fcTO**c4»n.  Hence, 
e(kel,/*,'*1'fm  for  ell  4  and  for 
ccmin(  (A/ ,  km*  ).  Set  mil  In  the  min  expres- 
•ion  to  obtain  a  value  for  e  that  is  independent  of  m . 

Tot  width  of  each  side  of  this  layout  after  the  intro¬ 
duction  of  drivers  is  tK+*n*n(m1/<,"*,)=n(n,'<,'»)).  and 
the  required  result  for  the  area  increase  follows.* 

Next,  we  consider  the  case  when  the  driver  is 
embedded  in  a  rectangle  that  is  not  square.  In  this  case, 
we  show  below  that,  in  the  worst  case,  the  longer  side  of 
the  driver  dominates  the  summation  in  equation  (2)  by 
proving  the  following  theorem: 

Theorem  3:  Assume  that  the  driver  for  wire  tu,  of  length 
l,  can  be  embedded  in  a  rectangle  that  has  length 
Oi=k,l{  along  tbe  direction  of  the  wire  and  width  bi=ktl{ 
ir.  the  perpendicular  direction.  If  0<s,f<l,  then  the 
worst-case  increase  in  the  area  A  of  a  VLSI  layout  on  a 
square  grid  is  max(8(A1/',,■*,).0(A1/<,■,*))  when  such 
drivers  are  introduced  at  the  beads  of  long  wires. 

Proof  Since  the  driver  is  no  longer  symmetrical  in  the 
twu  directions,  we  must  now  modify  the  equations  for  tbe 
ft  * 


«.=/  i.(«i . “m-*i . bm)=k,(L,*^(e,vai*i,^b/))‘ 

*>,=/»< a, . an.b, . bm)*kt(l,*  ♦«*,*>,))'.  (3) 

where  the  domains  of  /  It  and  /„  ore  extended  to  R*“  es 
in  the  previous  case.  The  d,y's  and  the  d^  's  are  deter- 
m.ned  as  before  (see  figure  9).  except  that  the  first  sub¬ 
script  of  1  refers  to  lines  parallel  to  wire  t.  and  that  of  2 
refers  to  lines  perpendicular  to  it.  Thus,  in  figure  9,  we 
have  the  following  values:  dlsjsl.  4a^x0.  4|*,=0.  4HJ*1. 
i„is  1.  ijcj  =  0,  and  the  other  values  are  all  xeros.  We 
now  require  tbe  fixed-point  of  the  function 

f(a  I . “m  • b  I . Am  )*(/  U.f  11 . /  Im  ./tl./ff.-  •••/ Cm  )' 

(where  the  arguments  on  each  on  tbe  right  hand  side 

are  ol,  .  .  ,in.b| . b„)  Again,  if  each  long  wire  is 

B(n')  long  to  start  with,  we  obtain  from  Corollary  1.2, 
the  following  upper-bound  for  the  width  of  each  driver: 

/?*max(e(n"/:,->).e<n’*'<'-«>)) 

■©(n**'  where  p  >mai(i  ,t ). 

As  before,  the  worst  case  occurs  when  r*l,  and  this 
gives  us  a  maximum  width  of  0(n*/<,"*))  for  each  driver 
along  its  longer  side.  We  still  may  need  0{n)  cuts  of  this 
width  in  each  of  the  two  directions  to  introduce  tbe 
drivers,  and  this  gives  us  the  required  result  for  the 
maximum  increase  in  area.* 

This  bound  is  once  again  tight,  and  tbe  lower  bound 
can  be  proved  using  the  layout  of  figure  3  with  n  /  2  verti¬ 
cal  long  lines  added  on  with  their  beads  along  a  diagonal. 
Note  that  in  this  case,  •  or  f  may  be  greater  than  1/2. 
and  so  the  area  blow-up  will  be  more  than  tbe  square  of 
the  original  layout.  In  particular,  if  p  a  1,  tbe  above  bound 
ceases  to  hold.  (This  will  be  tbe  ease  for  a  long  driver 
whose  area  is  proportional  to  the  length  of  the  wire  it 


drives,  i.e..  the  natural  design  for  a  standard  driver.)  In 
this  case,  the  convergence  of  tbe  system  of  equations  (3) 
for  specific  layouts  depends  on  the  constant  A,.  If  A,bl. 
then,  clearly,  drivers  cannot  be  introduced  in  any  layout 
containing  a  long  bidirectional  wire.  In  fact,  for  any  fixed 
value  for  A,,  we  can  design  layouts  in  which  drivers  can¬ 
not  be  introduced,  as  we  show  in  the  following  theorem. 

lbeorom  4:  If  tbe  driver  of  a  wire  of  length  (  must  be 
embedded  in  a  rectangle  whose  length  is  kl.  for  some 
constant  k,  then  there  exist  layouts  in  which  drivers  can¬ 
not  be  introduced  on  all  long  wires. 
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figure  10 

A  layout  to  prove  Theorem  4 

Proof  Consider  the  graph  of  figure  10.  It  is  the  same  as 
the  one  in  figure  3.  except  that  all  long  wires  are  b.dtrec- 
tional.  We  will  prove  that,  for  nr 2 /A,  drivers  cannot  be 
introduced  on  all  long  lines  in  this  layout  As  before,  let 
ua  write  out  the  equation  for  eech  d,  Tor  instance. 

di*k(l|  +  dm. ,♦<*». ♦<*») 

*»**(1i*rfi 


(where  m*(n/ 2)),  and  the  equations  for  tbe  other  d,'s 
can  be  written  out  by  inspection  Since  this  is  a  system 
of  linear  aquations,  this  can  be  rewritten  in  matrix  form 

as 

Ad -kl 

where  A  is  an  nxn  matrix  of  the  coefficient*  of  the  d('s,  l 
is  a  vector  of  the  Vi.  »nd  d  is  a  vect  he  d,'*  We 

need  e  solution  of  this  system  of  equati.  .he  d, 

nonnegative.  For  this,  we  note  that  the  matrix  /  -n 

If -matrix  when  n<(2/A),  and  from  the  theory  o'  » 
matrices,  it  is  simple  to  prove  that  tbe  system  o: 
tions  has  a  solution  with  til  components  of  d  nonneg„..e 
only  when  n  <  2/  k  [BP79].* 

For  particular  layouts,  we  mey  still  obtain  conver¬ 
gence,  but  with  a  large  area  penalty  For  example,  the 
layout  of  section  3  (which  has  ns  mutual  interaction 
between  drivers  and  wires)  will  require  ares  0((lvfc)m) 
after  the  introduction  of  eucb  drivers  This  result  is 
easily  verified  from  tbe  recurrence  reletion  for  s,  which, 
in  this  case,  is  *,■(!+* )a,.i. 

We  pointed  out  in  section  4  that  the  driver  has  a  last 
stage  with  a  vary  wide  channel  that  dominates  the  ares 
of  the  driver  Thus  such  e  driver  will,  indeed,  be  very 
long,  ond  have  constant  width  The  introduction  of  such 
a  driver  in  Ihe  circuit  can  be  disastrous  in  terms  of  ares 
blow-up  (if  not  impossible)  However,  it  should  be  possi- 
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ble  to  redesign  the  channel  to  occupy  a  more  square 
region,  possibly  with  some  penalty  in  the  total  area 
required.  Thus,  if  the  driver  originally  required  area 
9(1*).  then  this  design  might  result  in  a  driver  of  length 
©{!•)  and  width  9(1')  with  9(!,*')>9(1*).  but  with 
0<s.t<l.  Ir.  such  a  case,  drivers  can  a livays  be  intro¬ 
duced  in  any  layout  of  long  wires  with  a  polynomial  blow¬ 
up  in  area  given  by  Theorem  3  (the  degree  of  this  polyno¬ 
mial  will  be  very  high  if  s  or  f  is  close  to  1).  In  practice, 
we  expect  driver  designs  in  square  areas  with  no  extra 
area  penalty,  and  this  results  in  the  squaring  of  the  area 
in  the  worst  case. 

Finally,  we  look  into  the  case  of  a  rectangular  VLSI 
layout  in  an  nxn’  grid  (n»n  ).  As  before,  if  each  long 
wire  is  1  units  long,  then  the  width  of  any  driver  is 
by  Corollary  1.8.  By  observation  1.  the 
number  of  wires  in  the  layout  is  0(nn'/l)w0(A/l'). 
where  A  is  the  area  of  the  layout  before  the  introduction 
of  drivers  If  l*0(n' ),  then  we  may  have  to  make  2n 
cuts  along  the  longer  side,  increasing  its  length  to 
0(n  js/U-a)),  and  similarly,  8n  cuts  aiong  the  shorter 
side,  increasing  its  length  to  0(n'  Hence  the 

area  increases  fron  nn  to  £7(nn'  (*»-'(>-J>))i:0(.4,/<'"*>). 
If  I  =  fi(n  ).  then,  since  the  number  of  wires  in  the  layout 
is  0(nn  /l).  the  number  of  cuts  along  the  longer  side  is 
also  no  greater  than  Ofnn'/I),  so  that  the  length  of  that 
side  increases  to  0((nn /l)!*'!'"*').  The  increase  in 
length  of  the  shorter  side  remains  0(n'  l>,/<‘"*!).  Hence 
the  area  increases  to  no  more  than 

0{(nn  /  l)  n  0(A,'t,~p)),  since  l“0(n)  by 

observation  4.  Hence  the  worst  case  area  expansion  for 
rectangular  layouts  due  to  the  introduction  of  drivers  is 
no  grsater  than  that  for  square  ones 

To  summarise,  we  have  ahown  that 

1)  if  the  drivers  described  in  section  4  can  be  embedded 
in  e  square  reg.on  without  extra  area  penalty,  then  in  the 
worst  case,  the  area  of  a  VIS]  layout  grows  quadratically 
with  the  introduction  of  drivers  for  long  wires  This  worst 
case  occurs  when  drivers  require  area  proportional  to 
tht  length  of  the  wire  they  drive  (i.e.,  the  alandard 
design);  when  the  layout  is  in  an  nxn  square  and  has 
8(n)  long  wires  of  length  8(n)  each;  and  when  the  intro¬ 
duction  of  drivers  requires  9(n )  cuts  of  width  8(n  )  in  the 
layout  in  each  of  the  two  directions. 

8)  if  a  driver  requires  a  rectangular  embedding,  then  the 
worst  case  increase  in  area  is  8(.4I/'1-*'),  where  the 
longer  tide  of  a  driver  rectangle  for  a  wire  of  length  1  is 
8(1*)  grid  units  long  Ifpkl,  then  there  exist  layouts  for 
which  drivers  cannot  be  introduced 

7.  Application  to  Other  Layout  Problems 

The  technique  developed  in  the  previous  two  sec¬ 
tions  can  be  used  in  other  related  problems  in  VLSI.  Con¬ 
sider  the  following  question: 

Given  a  VLSI  layout  i uhal  it  (At  global  tjjtcl  of  making 
many  local  fpw.  »'  expansions  tAal  mutually  interact'’ 

The  Introduction  of  drivers  for  long  lines  is  one 
example  of  this  problem.  Another  example  is  fault- 
tclerant  computing  in  VLSI.  Von  Naumann[vNS6]  bat 
suggested  a  method  of  constructing  reliable  combine- 
tonal  circuitry  using  components  that  have  a  certain 
probability  of  failure.  The  construction  requires  replica¬ 
tion  of  tbe  basic  components  many  times. This  method 


does  not  carry  over  directly  to  VLSI  circuits,  since  the 
connecting  wires  here  are  also  likely  to  fail,  and  the  pro¬ 
bability  of  failure  increases  with  the  length  of  the  wire 
This  probability  can  be  reduced  by  increasing  tbe  width 
It  is  generally  accepted  that  a  fault  can  occur  aryvbere 
on  the  surface  of  a  silicon  wafer  with  a  constant  probabil¬ 
ity,  and  thus,  the  occurrence  of  a  fault  can  be  modeled 
as  a  random  variable  with  a  Poisson  distribution  Vr.der 
this  model,  it  can  be  ahown  that  the  probability  of  failure 
along  a  wire  can  be  bounded  by  any  constant  t>0  by 
making  the  width  of  the  wire  proportional  to  the  log  of 
the  length.  But  if  the  width  is  increased  after  the  circuit 
is  laid  out  (and  this  is  tbe  only  logical  way.  since  the  the 
lengths  of  wires  in  the  layout  cannot  be  known  before  the 
circuit  is  laid  out),  then  this  will  result  in  increasing  the 
lengths  of  other  wires,  and  it  is  not  clear  what  the  global 
effect  of  many  such  local  transformations  is 

We  analyze  the  general  problem  by  assuming  that 
in  the  worst  case,  a  "cut”  (as  described  in  section  2) 
would  be  required  to  make  any  expansion  of  one  dimen¬ 
sion  of  an  element  in  the  layout.  Thus,  to  introduce  a 
driver  of  length  l  and  breadth  b  at  the  head  of  a  wire,  we 
would  require  a  cut  of  width  f  perpendicular  to  the  wire, 
and  a  cut  of  width  b.  parallel  to  tbe  wire,  both  at  the 
head  of  tbe  wire,  in  order  to  accommodate  the  driver  Of 
course,  in  particular  layouts,  it  may  so  happen  that  com¬ 
ponents  are  sparse  at  this  region,  and  so  the  driver  car, 
be  introduced  as  it  is  without  any  expansion  of  the  lay¬ 
out  But  we  are  looking  at  worst-case  situations  and  »  e 
wii)  thus  assume  that  the  introduction  or  any  element 
can  be  achieved  only  by  making  cuts  of  the  appropriate 
widths.  Similarly,  for  wide  wires  in  fault-tolerant  com¬ 
puting.  the  width  of  a  wire  is  increased  to  u  by  making  a 
cut  of  width  w  parallel  to  the  wire. 

Our  general  construction  is  as  follows 
0)  Lay  out  tbe  circuit  (or  assume  that  the  circuit  layout 
is  given) 

1)  Number  the  spatial  elements  thel  are  to  be  varied 
(usually  this  is  s  subset  of  the  lengths  and  widths  ol  the 
lines,  or  the  dimensions  of  processing  elements  or 
drivers,  but  could  conceivably  include  other  objects)  as 
v,.  vj.  ,vr  in  some  arbitrary  order. 

8)  For  each  spatial  element  v,,  identify  tbe  other  ele¬ 
ments  Vj  that  are  made  to  expend  by  its  presence  For 
each  such  element,  set  6^-k.  where  v;  is  made  to 
expand  k  times  by  the  presence  of  v, 

3)  Tbe  following  system  of  equations  defines  the  final 
values  of  tbe  vt's: 

where  l,  is  a  constant  corresponding  to  the  imtia1  value 
of  some  linear  dimension  in  the  layout  and  /,  deter¬ 
mines  the  functional  dependence  of  v,  on  the  v, 's  If  for 
•acb  tbe  behavior  function  p,  is  sublinesr  with  break 
point  r,.  then  an  upper  bound  for  tht  v,'s  is  given  b>  R 
where 

R* max  r,. 

Examples: 

1)  For  drivers  embedded  in  a  square,  v,  is  the  length  of 
•  side  of  the  driver  of  the  ith  line,  and  i,  is  the  length  of 
tbe  <th  line,  and  /((x)sfc  **,  for  some  constant  k  and 
(or  aome  p  <  1 
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2)  Tor  fault-tolerant  computing,  v,  is  the  width  of  the 
ith  line,  it  is  its  length,  and  /,(x  )«A  In  *.  for  some  con¬ 
stant  k. 

B.  Conclusion 

In  this  paper  we  have  analyzed  the  effect  of  intro¬ 
ducing  drivers  to  speed  up  signals  along  many  long  wires 
in  a  VLSI  layout.  We  have  shown  that,  under  all  but  the 
most  naive  of  designs  for  a  driver  (i.e.,  the  case  when 
drivers  have  constant  width,  and  occupy  area  propor¬ 
tional  to  the  length  of  the  wire  they  drive),  these  drivers 
can  be  introduced  with  only  a  polynomial  increase  in 
area.  With  the  additional  assumption  that  drivers  can  be 
embedded  in  a  square  region,  we  have  shown  that  the 
area  at  most  squares  by  their  introduction.  All  results 
have  matching  upper-  and  lower-bounds.  We  have  also 
shown  an  area-delay  trade-off  in  the  design  of  drivers, 
and  we  have  generalized  the  upper-bound  proof  tech¬ 
nique.  Some  open  problems  that  remain  are: 

1)  Given  a  VLSI  circuit,  is  it  possible  to  design  a  layout 
for  it  so  that  the  presence  of  many  long  wires  is  "minim¬ 
ized"  according  to  a  suitably  defined  criterion?  Some 
work  has  been  done  on  minimizing  the  length  of  the  long¬ 
est  wire  in  a  layout  for  certain  classes  of  graphs.  [PRSBl] 
and  [RSSl]  have  derived  bounds  for  trees  and  [BL82] 
have  extended  the  result  to  classes  of  graphs  with 
separator  n,/s'*,c>0  or  Vn  log*n,tt0.  Leiserson's  lay¬ 
out  technique  [LBO]  minimizes  the  maximum  edgelength 
for  classes  of  graphs  with  separator  «'.  r>l/2.  Some 
lower-bound  results  are  presented  in  Leighton  [LBl], 
However,  our  results  on  drivers  indicate  that  it  is  not  the 
presence  of  a  single  long  wire  that  causes  the  worst-case 
increase  in  area,  but  rather  the  presence  of  many  long 
wires  in  e  configuration  that  requires  many  cuts.  Even 
heuristics  to  prevent  such  undesirable  configurations  of 
long  lines  should  be  very  useful. 

2)  Is  it  possible  to  incorporate  the  presence  of  drivers  in 
our  layout  model,  for  example,  by  using  nodes  of  variable 
size,  and  can  we  obtain  area-efficient  layouts  under  this 
model?  Using  such  a  model,  we  can  re-examine  the 
upper-bounds  that  have  been  obtained  for  the  layout 
area  of  many  common  circuits.  In  particular,  upper- 
bounds  for  AT°.a>0 ,  should  be  re-examined  under  the 
context  of  area-expansion  caused  by  drivers  of  long 
wires. 

3)  The  numerical  solution  of  equation  2:  In  section  6,  we 
used  the  Brouwer  fixed-point  theorem  to  prove  results 
on  the  existence  and  worst-case  behavior  of  the  solution 
of  equation  2  under  various  assumptions.  However,  the 
actual  solution  of  the  system  of  equations  for  particular 
layouts  still  needs  to  be  investigated. 

4)  Modeling  poly  lines  in  MOS  technology:  We  have 
modeled  the  wire  as  a  purely  capacitive  load.  However, 
poly  lines  have  rather  high  resistance,  and  may  be 
better  modeled  as  distributed  RC-network*  At  present, 
the  best  that  can  be  done  for  such  structures  is  to  use  a 
driver  of  constant  sire  after  each  fixed  Interval  of  the 
wire  and  this  reduces  the  delay  from  9(1*)  to  0(1),  where 
I  Is  the  length  of  the  poly  line.  It  erill  be  very  useful  to 
find  a  better  way  to  speed  up  signals  along  such  wires. 
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